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This is the final report describing the research performed under the NASA 
Grant number NAG-1 -477 at Hampton University during the years 1984-92. 
During the grant period, four M.S theses, NASA technical reports, and 
papers in refereed journals were published. Most of the research work 
performed under this grant were concerned with strong interaction 
processes ranging from kaon-nucleon interaction to proton-nucleus 
scattering calculations. Research performed under this grant can in 
general be categorized in ;o three groups: (1) Parametrization of 
fundamental interactions (2) Development of formal theory, and (3) 
calculations based upon the first two. 

Parametrizations of certain fundamental interactions, such as kaon- 
nucleon interaction, for example, were necessary because kaon-nucleon 
scattering amplitude was needed to perform kaon-nucleus scattering 
calculations. Of cause it was possible to calculate kaon-nucleon 
amplitudes from the first principle, but it was unnecessary for the 
purpose of the project. Similar work was also done for example for anti- 
protons and anti-nuclei. Formal developments to some extent were also 
pursued so that consistent calculations can be done. 

Note that all research projects supported under this grant NAG-1 -477 
were done with the aim that the results and knowledge obtained from 
these projects will be used in the shielding calculations of long space 
missions. 



The following is the list of awards, papers in refereed journals and NASA 
technical memos/papers published during the duration of the grant, NAG- 1-477. 
The rest of this report consists of a collection of these papers. Since the NASA 
technical papers, memos and students* M.S theses are very lengthy, we include 
only the cover pages of these. 
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Systems Division Best Paper for 1990-1991." 

Outstanding Publication within the Space Directorate at NASA/LaRC (1992) 
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“2nd Annual HUGS at CEBAF Proceedings", CEBAF Publication (1987), W.W. 
Buck, Editor 

"BRYNTRN: A Baryon Transport Computer Code: Computational Procedures 
and Data Base" J.W. Wilson, L.W. Townsend, S.Y. Chun, W.W. Buck, F. Khan, 

F. Cudnotta, NASA TM 4037 (June 1988) 

"Nucleon Nucleus Interaction Data Base- Total Nuclear anil Absorption Cross 
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"Relativistic Antinucleon-Nucleon Interaction: A New Level Ordering", W.W. 
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“The Third Annual HUGS at CEBAF Proceedings" CEBAF publication, W. W. 
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"Eikonal Solutions to Optical Model Coupled-Channel Equations", F. Cucinotta, 
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"Nuclear-Fragmentation Studies for Microelectronic Applications", D.M.Ngo, 

J. W.Wilson, W.W.Buck, T.N.Fogarty, NASA Technical Memorandum #4 143(Nov. 
1989) 
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" Meson Exchange and the Relativistic Multiple Scattering Fromalism", F. Cross 
and K. M Maung. Phys. Lett. Vol 229 188(1989). 

"Pseudoscalar pi-N coupling and relativistic proton-nucleus scattering", F. 
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"Progress in Proton Transport Code Development: Microelectronic Application", 
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"Kaon-Nudeus Scattering Calculations", B. Hong, W. W. Buck, K. M. Maung, J. 
W. Wilson, and L. W. Townsend, pl46, NASA/HBCU Space Science and 
Engineering Research Forum Proceedings, March 22-23, 1989, Y, D. Saunders, 
Y. B. Freeman, M. C. George Editors. 

“Covariant Multiple Scattering Series for Elastic Projectile-target Scattering" 

K. M. Maung and F. Gross. Phys. Rev. C Vol 42, 1681(1990). 

"Relativistic proton-nucleus scattering and one-boson exchange models", K. M . 
Maung, F. Gross, J. A. Tjon, L. W. Townsend and S. J. Wallace. Phys. Rev. C 
VOl 43 1378(1991). 

“Theoretical angular distributions from coherent subthreshokl pion production" 
P. A Deutchman and K. M. Maung. Phys. Rev. C. Vol, 41 R423 (1990). 

"A Nuclear Fragmentation Energy Deposition Model", Due Ngo, J. W. Wilson, T. 
N. Fogarty, W. W. Buck, Transactions on Nuclear Science, Vol38, no. 1, February 
1991. 

"Confining Potential in Momentum Space", J. W. Norbuty, D.E. Kahana and K. 
M. Maung, Can. Jou. Phys. Vol 70, 86 (1992) 



Seminars and Conferenc e Talks under NAG-1-477: 


"Optical Potential Calculations of Antideuteron Absorptive Cross Sections", 
Presented by W.W. Buck.with J.W. Norbury, L.W. Townsend, J.M. Wilson, 
Fall 1985 Meeting of the APS Division of Nuclear Physics, Pacific Grove, CA 

"Antimatter-New Frontier in Science”, Battelle Pacific Northwest Labs (March 
14, 1986), Invited 

"A Search for Galactic Anti-ion Cosmic Rays", 1st NASA-HBCU Forum, Atlanta, 
GA (April 20-22, 1986).Invited 

“In Search of Anti Iron", Florida International University, Miami, FLA ( 1987), 
Invited 

"Relativistic Isopin Exposure", The University of Kentucky, lexington, KY (May 
1988), Invited 

"Nuclear Physics at Hampton University", Virginia State University, Petersburg, 
VA (May 1988), Invited 

"Relativistic Antinucleon-Nucleon Interaction: A New Level Ordering", 
Workshop on Relativistic Nuclear Many-Body Physics, Columbus, OH (June 
1988), Invited 

"Antimatter Predictions from an Optical Model", The Ohio State University, 
Columbus, OH (August 1988),Invited 

"Relativistic Antinucleon-Nucleon Interaction: A New Level Ordering", Fall 
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Meeting, VA Beach, VA, Oct 21,1989 

Session Chairman, Fall Meeting Southeast Section of the APS, Tuscaloosa, ALA, 
Nov. 9, 1989 

"Nuclear Physics at Hampton University", Morehouse University, Atlanta, GA, 
Nov. 14, 1989, Invited 


"Physics: Science or Art", 1990 Stone Symposium on Sociology of 


Subjectivity. St. Petersburg Beach. FLA Jan. 25-28. 1990. Invited 

“Relativistic Dynamics: NNbar and qqbar". University of Illinois. Champaign. 
111. March 26. 1990. Invited 

"A Covariant Quark Model of the Pion". Argonne National Laboratory. 
Argonne.lll. Oct 16.1990. Invited 

“A Covariant Quark Model of the Pion". Carnegie- Mellon University 
Pittsburgh. PA.. March. 1991. Invited 
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A envarum fanmriatiiw nf I be nakii^c ital Kr ng serin far Ihr optical potential » prtvolol 
Wc consider Ike w nf a walir - uu ct c ou" MMcractmp with a vf>w zero twnpia zero .t-Wj target 
ibrtH ^h noon exchange. We On- - that a ctw jrium l equation lor the p r o pel i lc-targct r matrix can 
be obtained that ants I be ladder and dowsed ladder ihigrmr efficiently. From thrs equation. a 
nmhiptc scattering reran lor the optical p otential is derived, and *c show that in the impxJsc ap- 
proximation, the two-body f matrix r-svocialcil with the first-order optical potential is the one m 
winch one particle is kept on mass s.,dL The me a ning of r ati o ns terms in the mn h i p le scatt eri n g 
series is given and we descri b e bow to const met the first-order optical potential lor ciaslic scattering 
calculations. 


L INTRODUCTION 

It is well known that the relativistic (Dirac equation) 
calculations give superior results to the nnnrelalivisiic 
(NR) calculations in the case if p-nuclctts clastic scatter- 
ing and have been widely used in recent calculations. 1 
The first-order optical potential used in these calculations 
is the relativistic impulse approximation 1 (RIA) which is 
a relativistic ge n e r a lization of the noorclalivtslic //> ap- 
proximation. In the NR theory the ciaslic scattering of 
the projectile from the target nucleus is described by an 
effective interaction (optical potential) which is to he used 
in the Schriidingcr equation, and the scattering observ- 
ables arc then obtained. The optical potential itself can 
he expressed as an infinite series of scattering terms, sin- 
gle. double, etc., scatterings (hence the name multiple 
scattering scries) in which there are no two successive 
scatterings from the same target particle. Ity keeping 
only the first term of the infinite scries of the optical po- 
tential wc obtain the first-order optical potential. The tp 
approximation is achieved only after two more approxi- 
mations. namely the impulse approximation which treats 
the struck target nucleon as though it were free, and the 
factorization approximation which assumes that the 
range of the interaction is small compared to the size of 
the nucleus. The Iasi approximation is usually applied in 
order to avoid the complexities of performing the folding 
integral to obtain the optical potential. The existence of 
a multiple scattering senes for the optical potential (in 
fact there are several in the literature) provides tis with a 
means to calculate systematic corrections to the first- 
order results. 

In the relativistic p-nuclcus scattering calculations the 
effective onc-hody equal ion is taken to be the fixed energy 
Dirac equation. This choice is intuitively appealing as 
long as one considers the proton as an elementary fer- 
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mimi. hut the actual validity of this assumption is still ia 
question. This type of question will be answered only 
when one has the non perturbative aspects of QCD under 
the same de gr ee of control as in NR theories. Now the 
question arises as to what e ff e ctive intera c tion (optical 
potential) should he used in the Dirac equation to de- 
scribe /> -nucleus scattering. As mentioned above, in all 
the relativistic calculations, the optical potential used has 
been iiic RIA. 1 The RIA optical potential is obtained by 
simply folding a relativistic NN amplitude with the nu- 
clear density matrix. Strictly speaking, use of the RIA is 
an intuitive guess guided hy the nonrelativistic multiple 
scattering formalism, since a relativistic multiple scatter- 
ing theory (RMST) has not been available. 

It is important to realize that without a multiple 
scattering theory the f matrix associated with the first- 
order optical potential cannot he unambiguously deter- 
mined and consequently the characters nf the corrections 
to he made to the first-order optical potential are not well 
defined. The absence of such a theory prevents us from 
making systematic corrections, such as Pauli blocking, in 
a consistent manner. Therefore it is highly desirable to 
have an RMST. Probably the most appropriate approach 
might he to apply the methods of field theory to the prob- 
lem. But the development of the RMST in this direction 
has been hampered by the problems arising in the treat- 
ment of the interacting manv-hody ground stale,' 
description of I Ik- residual interaction he! ween the projec- 
tile panicle and the target constituents, and many other 
obstacles not encountered in the NR theory. 

In this work wc lake a less formal hut more intuitive 
approach and describe the projectile-nucleus scattering 
problem in a meson exchange model. A brief account of 
this work has already been given. 4 In this paper wc will 
develop the ideas reported there in more detail. 

l-'rom the beginning, wc would like to make it clear 
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aim is lo derive a multiple scattering theory for 
the projectile-nucleus scattering in the contest of meson 
exchange- We will not consider the full implication of 
the formal field theoretical treatment of the interacting 
many-body problem, which is admittedly very difficult . 
We develop an approach which permits the standard 
multiple scattering techniques of NR theory to he applied 
with slight modification. We ignore the full complica- 
tions of antisymmetry required by the Pauli principle and 
ahu ignore the spin pa.t of the problem. We do not as- 
sume any particular form of equation for .he projectile- 
nucleus r matrix, hut begin with the most obvious tact 
that the I matrix is obtained by summing an infinite set of 
diagrams in which the projectile is interacting with the 
target panicles through meson exchange. Although the 
approach we lake in this mirk seems much less compli- 
cated than the format field theoretical approach, it has its 
share of problems, such as the appearance of spurious 
singularities and the necessity for judicious treatment of 
the crossed meson diagrams. In this work we show Iriw 
an RMST can be formulated within the context of meson 
exchange, unambiguously determine the r matrix associ- 
ated with the optical potential, and show that under the 
impulse approximation the I matrix lo he used in the 
first-order optical potential is the solution of a rdalivrcic 
two-body equation in which one panicle is kept on its 
mass shell. 

This paper is arranged as follows. In Sec. II the deriva- 
tion of a multiple scattering series of the optical potential 
is reviewed, following the approach of Watson.' In Sec. 
Ill A it is shown how the crossed meson diagrams should 
be treated together with the box diagram for the case 
where the intermediate state target is in the ground slate. 
In See. Ill B the complications that arise in the case of in- 
termediate target excited stales arc discussed, and it is 
shown bow to handle the box and crossed box diagrams 
in this case. In Sec. Ill C a covariant equation for the 
projectile-target r inalrix is presented and it is shown how 
a multiple scattering scries for the optical potential can 
be derived. The meaning of various terms in llic multiple 
scattering series arc discussed, and it is shown that the 
most appropriate two-body r matrix to be used in the 
first -order optical potential under the impulse approxima- 
tion is the one calculated from a covariant equation in 
which one particle is kept on the mass shell. A general 
discussion and conclusion follows. 
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H. NONRKI-ATIVISnC FORMALISM 

In this section the nonrclaiivisiic multiple scattering 
formalism is reviewed. The approach of Watson' is fol- 
lowed, since it is more closely related to our relativistic 
formalism than the more commonly used Kcrmaii- 
McManus-Tlialcr* (KMTl method. Since we arc interest- 
ed in deriving a relativistic multiple scattering series for 
the optical potential, we will bypass the multiple scatter- 
ing treatment of the r matrix and will concentrate, in the 
present section, on the multiple scattering analysis’ of the 
nonrelalivisiic optical potential. 

We begin with a total Hamiltonian II for the 
projectile-nucleus system given by 


"'*2 

//=//.+ K. (2.1) 

with 

//„ = //,+*„ € 2-21 


ami 

"«= f *.+ !>., • OJ) 

.i .- / 

V = 2 r. - 12.4) 

* • 

Notice that the total Hamiltonian II is separated into two 
parts, the unperturbed Hamiltonian //„ and the residual 
interaction Y. It is the separability of II into ll u and K. 
which permits the derivation of a multiple scattering for- 
malism. The residual interaction Y rs taken to be the sum 
of the two-body interactions between the projectile parti- 
cle "IT and the target panicle ”i.“ The unperturbed 
Hamiltonian //„ is written as the sum of the target Ham- 
iltonian II , and h a . the kinetic energy operator for the 
projectile. In NR formalism, the target Hamiltonian is 
just the kinetic energy operators of the target panicles 
plus the sum of their pair interactions. 

The separation of the total Hamiltonian in Eq. 12.1) im- 
plies that we have some means in finding the solution to 
the target Hamiltonian II Therefore in NR theory the 
complexities of the d-budy problem arc separated from 
the rest at the very beginning. Now write the 
Lippmann-Schwinger equation for the projectile-nucleus 
r matrix in operator form as 

T = V + VG a T 12.5) 

with 

G„ - 1 E -If.+ iij) 1 . (2.6) 

Here G u is the unperturbed Green's function and the i »| 
prescription has been used to incorporate the outgoing 
boundary condition. The many-body nature of Eq. (2.5) 
is apparent since the propagator G u involves the target 
Hamiltonian II , - 

I-or clastic scattering problems it is useful lo introduce 
an effective one-body potential (optical potential). The 
optical potential is defined as the potential that describes 
the passage of a projectile through the nucleus with the 
nucleus treated as a passive medium, i.c.. the nucleus is 
treated as though it cannot be excited, lo accomplish 
this, first detone a projection operator P which projects 
onto the ground slate of the target and Q which projects 
onto the excited slates of the target including the breakup 
slates. Therefore, 

r + Q - I , (2.71 

where 

= <2.«> 

and |<A„) is the target ground state. Now Eq. (2.5) can he 
rewritten as 

7 U ) UPG„PT . 12 •») 
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U - I' i l’QG a QU . 12.101 

The U appearing in these equations is ihe optical poten- 
tial operator and Eq. (2.41 together with Eq. 12.101 are 
equivalent In Eq. (2.51. 

Sincr we are dealing with strong interactions, it is im- 
practical to solve Eq. 12.101 for V as it stands. It is at this 
point that the multiple scattering approach provides us 
with a Mg advantage. We may express U as X.' i U . and 
rewrite Eq. (2.101 as 

V. -t, *r,QG„Q 2 U, (2.11) 

i i 

Now define the Watson r operator as 

r , ~ r . ' r ,QG a Q r , 

and observe that Eq. 12.1 II can he written in terms of r. 

U. = r.*r.QG m Q?.t', «•»» 

I ■ 

Summing over the index i in the last equation gives the 
Watson multiple scallering series for the optical potential 
operator. 

«/- i.r.QC m Qj.V, (2.14) 

f I * I #4 

Notice that Eq. (2.141 is an infinite series in r instead of 
the two-hodv interactions r as in lug 12.101. Each term of 
Eq (2.141 can he interpreted as single scattering, double 
scattering, and so on. hence the name multiple scattering 
series. Ify keeping only the first term of the series we ob- 
tain the first-order Watson optical potential 

(/•"-- Jr, . 12.15) 

• i 

The operator r is not the free two-hodv r matrix be- 
cause of Ihe many-hody propagator in Eq. 12.121. but re- 
lated to it by 

r I glr . t2.lf»l 

where iIk- free two-body /-matrix is defined as 

/ -' r I rgi , t2. 1 71 

wirh e the free two-body propagator. I -or high proieelile 
incident energies one usually approximates r by f (im- 
pulse approximation') and obtains for the first-order 
Watson impulse approximation optical potential 

i 

gri'* y ( . (2.IKI 

Im X* r 
r I 

We can also rewrite Eq. t2.l4t in terms of the free two- 
body / matrix, r. 

v £/. i £/,•<*;..(> riu. « £ LC’f.UN/, i2.i‘)i 

/ i • i 

As we have mentioned above Ihe first term in Eq. (2.1**) 
gives the fust -order impulse approximation optical polen- 
lial. The second term can be interpreted as Ihe propaga- 
tor correction term. Ibis term origiiiaics from the fact 


that we have written the optical potential in terms of the 
free f matrix r instead of r. For high projectile energies 
Ihe differences between / and r become negligible and the 
impulse approximation should give good results. The last 
term represents the multiple scattering terms. For NR 
scattering calculations the r matrix appearing in Eq. 
(2. I*>) can he obtained from Eq. (2.17) by employing a 
choice of r. for example (he Reid potential, or by filling 
the iViV experimental data directly by nsing an appropri- 
ate functional form. After the choice for the t matrix is 
made, solving Eq. 12.17: together with Eq. 12-91 is just a 
technicality. 

III. RELATIVISTIC FORMAIJSM 

In the last section we reviewed the nonrefcilivtsiie mul- 
tiple scattering formalism and outlined bow a multiple 
scattering scries for the optical potential can he obtained. 
We pointed out that Ihe key feature that enables us to 
construct a multiple scattering scries is the s epa r ab ility of 
the total Hamiltonian into an unperturbed llamilloman 
dcscriMng the free project ilc-largcl system and Ihe resid- 
ual interaction which is the sum of the Iwo-hndy interac- 
tions betw e e n the projectile and the target particles. Un- 
fortunately. there is no analogous procedure in the rda- 
livisli." case. Eirsl of all. one cannot naively write the tar- 
get Hamiltonian as the sum of the Dirac Hamiltonians 
plus the' sum of iwo-body interactions, since the Hamil- 
tonian written in this manner docs not have a lower 
hound. 1 In order to treat the projectile-target scattering 
consistently in a relativistic formalism one needs to resort 
to a field theoretical approach. 

In this work we lake a less ambitious route and show 
that j reialivislic multiple scattering series can be formu- 
lated in the context of a relativistic meson exchange mod- 
el. In the following wc will consider a scalar “nucleon” 
interacting with an .4 -body spin zero isoxpin zero target 
where the interaction between the projectile and the tar- 
get is described by meson exchange. Since we do not as- 
sume any particular form of equation for the projectile 
target / matrix, wc will stall from the most obvious fact 
that it can he obtained by summing all possible meson ex- 
change diagrams for I Ik- projectile large! system. A 
minimal set of me son exchange diagrams required for any 
such theory is the set of ladder and crossed ladder dia- 
grams. In ilie limit when the heavy target becomes 
infinitely massive. Ihis set reduces to a one-body equation 
for Ihe lighter particle moving in an instantaneous polcn- 
lial produced by Ihe heavier particle tlhc onc-hody lim- 
11 * 1 . and at high energies gives the ciknnal approximation 
to scallering.'' In this work wc seek a theory in w-hich 
these relativistic ladder anil crossed ladder diagrams arc 
summed efficiently. 

In Eig. I l lie target is represented hy a double hue. the 
dashed lines represent the exchange particle < meson!, anil 
the solid line repieseiits ihe projectile. For the intermedi- 
ate states the largei can be m Us ground stale. denoted by 
n --(i. or in cxcilcil siaics. «••(>. which incluilcs ihe 
breakup stales. The nolaliou is very compact: each dia- 
gram in Fig I actually represents a set of diagrams which 
can Nr ohia<ncil bv opening up ihe bubbles ai the meson- 
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FIG. I. The pnyed ilMarjrt / noiriv n Jum (he son rf all wra ciduagc pruccaes up to Ike sit aider dugramv The m- 
jle Ink iqwnnte the projectile am] the doable hoc represent* the target. The dashed lines are the exchanged nr enm . (a) is the one 
meson exchange term. Ibl is the box. and lei the crossed box. In the fourth and higher order diagrams, all p oss i b l e intermediate target 
stales are sommed. 


target vertices. For example, the set of diagrams con- 
tained in the box. l-lg. I Ibl. and crossed-hox. Mg. Ilcl, are 
shown explicitly in Figs. 2 and J. 

In our view, the solution of the relativistic problem in 
the meson exchange approximation is equivalent to 
finding an integral equation which sums all of the dia- 
grams shown in Fig. I. The construction of such an 
equation confronts us with three problems. The first 
problem, which does not occur in the nonrclalivislic case, 
is the ap|iearancc of the crossed meson diagrams. These 
and all oilier irreducible diagrams (i.c.. those which can- 
not be separated into two pieces by a line which intersects 
only the projectile and the target) will he included in the 
kernel of the integral equation. The second problem con- 
cerns the I real men i of excited states. All diagrams, ex- 


cept for the one meson exchange term, include terms in 
which the target propagates in an excited slate. A third 
problem is that each diagram includes terms in which the 
projectile may interact with two or more different target 
particles i.nulltple scattering). In this section, we arid 
first discuss how ilk' crossed diagrams arc handled, and 
then discuss I lie complications arising from the oc- 
currence of ext led stales. 

A. C—rcHatMwi b et ween the box anti twnd-Wt d i x g rxu rt 

We know from the two-body problem that the ladder 
sum does not give a good approximation to the true solu- 
tion of the Bethe-Salpeicr equation. There is no reason to 
believe that it would he otherwise here. In fact, in the 



FIG. 2. Figure Kb) is redrawn by opening up l lie bubbles al I tie meson-target vertices. The sum is over the target particles. 



42 


OOVAKIAN1 MULTIPLE SCATTERING SERIES FOR ELASTIC . . 


1685 



l it •- 5. figure 11(1 is redrawn hy «*pcnitig lip I Ik hubbies al the moonaargei vcrlicw The Mtm IN oirr the larpl p-nielcN. 


Ivivhndv priWem I Ik- box diagram and the mi\MxM>n\ 
diagram tend la caned, showing llial it k uii|irslificd in 
ncgkxl crossed racyw diagrams. 

In this see Ik m we show llial the cancellations hdwcvn 
l he hnx diagram ami I Ik* cross-box diagram still occur in 
Ihc case of pri^vlik- nnclcns valtcring. In nnlcr in 
d c m nns lralc this canccllalinn. «v perform ihc integral inn 
over ihc relative energy for the intermediate stales. 


In Mg. 4 the hnx diagram. Fig. Itb). and the crosscd- 
hnx diagram. Fig. He), arc redrawn and for the fnar- 
mnmentum variables fnr each internal partK'lc are la* 
heled. l : nr nur purpose now, it is sufficient In consider 
the target as a structureless panicle as shown in Mg. 4. 
Ilv employing standard Feynman rules the hit and 
erossed-bnx diagrams nT Fig. 4 are given hy the following 
expressions: 


JW'"- 




d X k'dp'„\v*' - let k I p „ I - »i/J - 

~r» *»l ]( F:;»k'i -(»’ I - * »» 1 


.w 


tin* . 



/ 


d'k'dp'„\, ' lelkl - p'„) : ii|I - 

fc’tq) - (2elk> /»')-" /»/ HCjtk’l ' *»' P« it/] 


where tin: total four momenta in the e.m. is 


an 

13.21 


/. hp--i ir.oi P ’ » r -/»" » r" . 

and the three monK'nta and the on-shell energK*s are ilefineil as 
p=-P k: q p p I p". p’ P‘ k’ . 
w = f/r4lk' k»‘j' etk* tm’ » kV T.ikl t.V; » kV * . 


We assmiH* forward scattering, i.o.. k - k”. so dial the inesnn poles become ilonble poles. The external particles arc lak- 
en to he on tlieir mass sIk-II. 

Mgnres 5tal and 5lbl show the locations of ilie poles twlien ikl ami ik’l are small* for the box diagram and c rosso I- 
box iliagram. respectively. We will etaliiale tlie hn ami erossetl-be-x diagrams by using the resnlne lho»re«n. In iIk* fol- 
lowing expressions the superscript on 3/ distinguishes between die box and the erossol-liox diagrams, the suliscript is 
for the type «»f ptile under eonsideralKni. and the letters I' nipper half planet ami /. Ifor lower half plane) aie usol to re- 
mind hs how the citutciir is closed. For example 3/ 1 Ill means the negative energy proieetile pole 1 pi contri- 

bution from the fonnh-or«ler box diagram (4.4 i for n * t), ami the irdegralion contour is closed in the upper half plane. 

Fvaltialc the box diagram for n tl. Close die contour iii the upper half plane and piek up the target positive energy 
pole, double meson pole, anil the projectile negative pole: 

.»/ "if'.n m in i .i/'. 'tf-.u n» i «' •*#?.« m . 

* r 


where uv have iimxI the subscript - suggest: wit for the meson double pole eontribulion. These •-oiiiiihulnms are 
-* - (/;,.ik i f: (1 i , <n'| : 



Ml 
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, i — j 
i 

M'/W.n 

Ml 
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i2r> 1 
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where 

B =2«(l:* 0 (kl+«)(r I (k')-|r(k»-»)’l . 


For I hr crosscd-hox diagram we close the contour in the lower half plane for all « and pick up the double meson 
poles and negative energy projectile and target poles: 

M** r C/....I+ . 

The individual pole contributions arc 




d'k’[F f G) 


d'k' 

I2r>’ |« -ie(k)Te(q»)-‘) > (2e(q)l|Aj(k't-(#: 0 (k)-e(k)-e(q))- , J 

4 

(2rl* 


M 4, AL.m)=-~ g — f — - 
11 12rl* J |w i 

.. K r klUf.Ck')) 1 ! 1 

M 4B ,U..*)= - - I ~T I — — . 

[e J (ql — -F„(kl-li.lk'l)-J|2/;.(k'II 


where 

F=l£;(k')-C£ 0 (k»-») ; J|[e I (ql-(e(kJ-«)-’H 2«u(eCkl-wM . 
G=2*XF 0 (k|-«>le : lq»-(e(k)-») ; | . 


At this stage one could show that, at threshold, the 
dominant contribution of the bos diagram for a — 0 
comes from the positive energy target pole and the meson 
poles give the second largest contribution. For the 


9 9 ' P" 


Jg, p, 

a 

• 

V 9 - 9 - 
• 

V P ’-P” 

r : p = p*p-p- 

p*p-p~ 
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\P''P" 

. 
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p A x 
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p+p-p“ 

X) 

■Q> 


(b) 

FIG. 4. Figures Kb) and ltd are redrawn with explicit labels 
for the projectile and the target momenta. 


I — — 

crossed -box diagram, the treson pole contribution is the 
dominant one and is nearly equal to the meson pole con- 
tribution from the box diagram but with a relative nega- 
tive sign. Since we arc interested in other energies beside 



p^ plmtxv (Croxsc4 Hul) 
2c «’ r x 




(t>) 

FIG. V Tlie pole structures of the box diagram [Fig. 4(a)) 
and the crosscd-box diagram |l r tg. 4(b)) are shown in the com- 
plex plane. 1 lie circled dots represent the double meson 
poles. 



threshold, we evaluate the various pule contributions td the leading terra, an order of magnitude htfrr Ihatt 
without any further approximations. I lie only restriction when the target n ow-shcll. 

is forward scattering. l-'tguie 7 shn»v the A dependence id these caneelta- 

Ftgure b demonstrates the eancellalton between the lionv The legend of the turves raean the same as tn 


line in it M * ' F t M\ * \/M *, *, I, tlie absolute magnitude id 
the ratio of the sum of negative energy projectile pole and 
meson pile eonlribution to the positive energy target pile 
coni nhut ion for n - 0. The dolled line is the ratio of all 
the pole conlrthiiiions from the crossed l*ox «•» the posi- 
tive energy target pole of the bos diagram, fhese two 
liras. lie praetieallv* on lop of each other. Finally, the 
solid lira* is ft II* J, t If** I If * *, 1 which is the ra- 

tio d the sum id the fufl erosseil bo* and negative energy 
projectile pole pins the meson pAs id the bos diagram to 
the positive energy pole ctf the bos diagram. In these cal- 
culations. the target mass is taken to he Af„ Hun. where 
m is the mass rd the projectile particle and the meson 
mass is taken to he ft ~m /l. Ibis figure shows (hat. 
when the target is in the ground state In - 01. tin* poles id 
the ho* diagram, winch remain after the target h, pirt on- 
sMl, and the erossed-hos diagram, are era'll of tin- outer 
td W>*- id tin* lead i ng M*,* t term, ami hence arc far 
from being negligible. However, when the hoc and 
crossed bos are taken together, an excel ten I cancellation 
occurs, as shown by the solid line fiwr the energy range 
shown. After the cancellation, the positive energy target 
pole clearly dominates, ami whatever is left over is less 
t han B.J'i of this dominant eontrihulion. 

If the projectile is put on mass shell, instead of (lie tar- 
get, the ratio of the correction from the box and the 
erowod hoc to the leading terra would he 
It A# 4 * » M* * I Jt#**1/Af*4L amt ihK ts the dra-dashed 
line shown in Fig. <*. This result shows that the cancella- 
tion between the box and crossed-box diagrams is not as 
einnplete when I lie projectile is mi-sht-lF Inti still ijuite 
grant. The terms winch remain are now Mravii I 4 r » 

n o 

I • ■ ! *■*•«■* ; *. ‘ J-" -■ § 



;iip» »>.•** T *>#*r» u»*f* 

% % m<~m 

I ll*. **. I hr rs««ri:ilv*Nt»tix fh'mwi I lie tw«* ilMtJMm diitl llir 
mx vh**w*t hn tt li, IIh Htfvi K 

t**km f** % Uw* k licit ’ m ** stir m^ss t *i iIh |w**tktf*l%\ Vr 

lt*e m iIh* txsi. 


,W„ Am. where the binding energy is neglected. The 
projectile laboratory kinetic energy is ltd at I GcV. As 
can he seen from the solid lira*, d the target is on-shell the 
cancellations become better as A increases, and exact 
cancellation occurs when A »r. A\ shown hy the solid 
line. Ihts is an excellent approximation even ft*r light nu- 
clei. If the projectile is on -shell the cancellation does not 
improve as A » r . , shown by tl*c dot-dashed line, 
reflecting the fact that, in this case. *hc cancellation de- 
pends on the properties of the projectile and ran era the 
target. 

Tin.* above results suggest that, when ihc target ts in 
the ground slate, it is an cudlnri approximation to keep 
only the positive energy target pole for the intermediate 
slates, which is equivalent It* keeping the ground state 
target on its mass shell in att intermediate stales. The 
cancellation is less complete and I lie approximation levs 
accurate for realistic eases with spin and charge ex* 
change, so that it is desirable, in the general ease, to in- 
clude tat lead in principle* these extra terms as higher or- 
der corrections to the kernel [they bertmte part of F* in 
Eq. t * 51 as dscrihnl below). 

The alternative approach td put ting tin- projectile on 
shell has been seen to be less rail justified: the additional 
correction terms arc brger and do ran decrease as 
A « « . We believe that this analysis provides a satisfac- 
lory motivation for wrong a fixed energy Dirac equation, 
in which the projectile »s od-shell anil the target is on- 
shetl. to dcscrilie dastie nueienn-nuctews scattering. 

It. Trralnwnl nt tlw c veiled states 

In this suliwclion we will consider bow to treat the in- 
termediate target excited states. It would he tempting to 


Tf.tb t OV it t) 



:•'» a* m 

M* -m ■■ f m 


t It!. ?, I he taijset maw dcpi’uth «*•*■ **i the ranu-Halioi* o 

shown for « tl I hr protcvltk' hilxirariuv linrIK'oer-i is i,«k 
i » I 1 ' K' I t id Vv »h- ihw iinsi. » t tit the test 
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say that I he state apftraaimatioa that we have advocated 
in the case of a — 0 shuold work here also, ami that the 
excited state of the target should be pat iw mass shell. 
Bat for a S O, farther complications may arise be c a u se of 
the so-called dissolutkai singularities . m 

The dissohitioa singularities are spurious singularities 
which arise when a highly excited heavy target is pat on 
its mass shed. To see hour they come aimat. consider pal- 
ling the excited target on its mass shed in the expression 
for the box diagram, i.c.. 


I«iik'i-iir-p;i-*-„| 1 - 2wi 2Kjk-) ■ 
The projectile propagator in (3.1) can be factorized into 


Ie ^ ^k*)-p; ^ ^ , ®Ie(k*l+CB r -Jl.lk•))^ , 


X[elk'l— (IK —folk'll) 1 . 


In the last expression we see that there are two singulari- 
ties, one at IK=rlk')+£„tk') which is the usual clastic 
cat and the other one at IK=£ - (k*)— e(k') which is the 
dissolution singularity. This second singularity is spuri- 
ous because it docs nut occur when the diagram is calcu- 
lated exactly. (It can be shown that it is cancelled by a 
simitar singularity in the M* * p term.) When n =0, this 
singularity occurs at R'-E,-r, which is way below 
threshold and hence not of importance. However, when 
the intermediate state is highly excited In / 0), the singu- 
larity can move into the physical region and is a cause for 
concern. It has been an obstacle in developing an RMST. 

To see when this singularity becomes potentially 
dangerous, study the locations of the pities in the box dia- 
gram. Fig. Sfa). By approximating IK ~ M a + m (thresh- 
old) and taking |k*| to be small so that flk'|a« and 
C a (k’l«lf a , we see that the negative energy projectile 
pole and the positive energy target pole arc separated by 
an amount M a —M m +2 m in the upper half plane. As M, 
increases, the positive energy target pole moves towards 
the negative energy projectile pole and when the excita- 
tion energy of the target reaches 2m the pules touch and 

I 


a singularity arises. In this situatkm it is dearly not a 
good approximation to lake one of these poles and 
"neglect** the other. In addition to these spurious singu- 
larities in the projectile propagator, there are other spuri- 
ous singularities arising from the meson propagators 
when the excited state target is put oa its mass shdL For 
catarialional purposes these meson singularities are even 
worse than the ones from the projectile propagator siace 
they can arise fur rdal : vdy low excitation energies. At 
threshold they will appear when the excitation energy 
reaches the meson mass. 

The situation in the lower half plane [Fig. Sfa)) is 
difermt. As E n increases the negative energy target pole 
moves away from the positive energy projectile pole. To 
sec this explicitly we put the projectile on its mass shell in 
Etj. 13.1): 


, , . «p;-e(k'» 

(' inti- 


and the target propagator is now 

[A'ilk*)-! IK -#»;>) , = [fs B tk'l — (IK — etk*l)J 1 


X[#:Jk')+(IK-e(k*))| 1 

and exhibits no spurious singularities in the physical re- 
gion. It can easily be seen that meson propagators do not 
have any such singularities either. 

The above analysis suggests that, when we evaluate the 
expression 13.1) for a =0, we should close the contour in 
the upper half plane (to obtain the best approximation), 
but for ft '/ 0 we should dose the contour in the lower 
half plane to eliminate the problem of spurious singulari- 
ties. 

We now study the accuracy of this prescription by 
evaluating the box diagram for n / 0 hy closing the con- 
tour in the lower half plane. The contributions come 
from the positive energy projectile pole, double meson 
pole, and negative energy target pole: 


M*\L.n+ 0)=M\*(L,n/0>+Afi*(f-.» * \\L,n /0 ) . 

The contribution from these poles is 

d'k * 


U* A It a / — _ X f 

' (2*)' J [w J — (etkl—r(k’)) J J*’[2c(k‘)KA**lh')— (IK— e(k*)l*J * 

^t.t t£. n '/ Q) — )» f 5^. A. 

' ’ ( 2 «rl* J |w- , -(li„(k)+is: i ,(k*)) J ] I le-’tk , l-(IK+fc*.(k')) J J[2h*.(k*)l 

/ 0 )~ - -**- f — - ; ^'AKJ rT , 

12*)* J 4m'( .• (k'l lrlkl l <-. ) ) J ( K'i k' > ~ i K„i k ) - <■> )* ) ’ 


where 


C**[lt-(k , l-(ff 0 fk)-«) I l|[e-’lk')~(t*Ck) + r«i)-’)-2«a(e(k)+w)| , 
D ^Zid^lkl-wlIr-’lk'l-Celkl+wl-’j . 
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In Fig. 8 calculations of the * ^ 0 casts arc shown. The 
dashed line is the ratio of the snot of the meson poles 
(box) and target negative energy pole (box! to the positive 
energy projectile pole (box), U., M 4 / 

The dotted line shows the ratio of the crossed -box dia- 
gram to the positive energy pole of the box diagram. The 
solid line is which is the ra- 

tio of the stun of the fidl crossed box and the negative tar- 
get pole (box) phis : icson poles (box) to the positive ener- 
gy projectile pole. For this calculation the target mass 
and the meson mass are the same as the n — 0 case, and 
the excitation energy of the target is taken to be 
Am =m/IOO. As in Ibe • *0 case, the canceRattons be- 
tween the bn and the crox<d box still occur to a very 
large extent, although the cancdlalinn is not as good as 
in the previous case. It can be seen that, after the cancel- 
lation, the leftover terms in the energy range shown are 
less than 4% of the dominant projectile positive energy 
pole rod Inbuilt xi. 

In Fig. 9 the curves mean the same as in Fig. 8. hot are 
shown as the function of target excitation energy with the 
projectile laboratory kinetic energy fixed at I GeV. As 
before the sohd fine shows the net result and it is seen 
that the cancellation becomes better as the excitation en- 
ergy increases. This is a good signature, since as the exci- 
tation energy becomes higher, the one particle knockout 
term, which is the dominant inelastic contribution at 
medium energies, will become more important, and it is 
advantageous to have the cancellations improve as these 
conthboiknis become larger. 

The preceding analysis suggests that for the intermedi- 
ate target excited states, keeping the positive energy pro- 
jectile pole of the box diagram provides us with a very 
good approximation and at the same time avoids the 
spurious singularities that would arise by putting the ex- 
cited target on the mass shell. 



IT* i. K. The cancel tot inns between tin- bos diagram ami the 
croxsed-hnx diagram are shown for n ' O. M„ is tin- same as 
Fig. h ami the excitation energy Sin - in /MXI. See the discus- 
sion in the text. 
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l-Ki. V. C'anerltolioars ftw • Ml shown as a functhm of esetia- 
lion energy. See the diseussiun in the text. 
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In the preceding subsections, we have discussed the 
cancellations between the box and the emsed-box dia- 
grams for both xr =0 and n / 6. We have concluded that, 
because of the excellent cancellations between these dia- 
grams. we should keep the target on the mass shell for 
« =0. and that for n * 0 we can avoid spurious singulari- 
ties and still have a very good approximation if the pro- 
jectilc is kept on mass shell. Note that these approxima- 
lions arc obtained by considering the N»x and crossed box 
together. If one considered the Nix diagram only, these 
approximations would mil lie as good, as can he seen by 
tire dashed curves of Figs. 6 and 8. 

We now follow the suggestion provided by the last two 
subsections and write an integral equation for Ibe 
project ilc-largct t matrix in the following form: 

7 = V +■ W,' pFT + F{K»; „Qf . |J.J| 

where (*/ „ is the propagator for the target in its ground 
state and on its mass sliell. fij „ is flic propagator For the 
target in its excited stales with the projcelile on mass 
shell, ami I* and Q are target grout’d state and excited 
slate projection operators, respectively. The three- 
dimensional propagators appearing in l-q. l.U) van Ik- 
written in a manifestly covariant form as 

’* " lie)' I nr p'’l 
H’lm’h in e.iw. frame is 




.’mil 


nr 


I *!/»;. (»' /:„lk'»J 

:E 0 tk‘>|. 'lk’) -III' -i; ft ik'i)’| * 

tt.Vi* 

) h ' Ip’’ til') 

{in)' \M- UF - p'l-’l 
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which in c.m. frame is 

I 




Mp' t , e ( k' I ) 


<2jrl‘ 2etk* »{/:,; *k*>~l IK «tk*»)‘| 


I.UId 


liquation t3 3| is I Ik* first major result of this paper. It 
is a three imei.sional covariani equation for llie 
project ile-iargct r matrix. If iIr- iiiuIimir-mhi rtrluMp* ir- 
reducihlc diagrams ito he described below! are ucglecvd 
the driving; term of this equation. I*. assuiiK-s a very sim- 
ple form. It is iIr- sum ofone meson exchanges between 
the projectile and the target particles. Hr- special feature 
of this equation is that it has two three-dimensional prup- 
Liptoni in wliR'h the target is .m mass sIr -11 when iIr* tar- 
get is in the ground state aod the projectile is on mass 
shell when the target is in an excited stale. 

lit llq. djl I' is iIr* sum o ' alt irreducible !«R*s>m ex- 
change contributions. In iwrilcr to derive a multiplc 
scallcring series I':h llie optical potential in a fashion 
similar to NR ihoiry we write 

l~' - VT-, if". t.V4l 

where the first term is the one meson exchange -liagrams 
summed over iIr* target panicle iiRlex. Hu* second term 
I" includes the contributions fnnn the irrediR'ible dia- 
grams involving multimcson exchange, To he explicit, 
the fourth-order contributions to I" would he the parts 
of the box diagram, which do in M have the target mi mass 
shell when u ~0. and the projectile on mass shell when 
it /• U, and tin* emssed-hox diagram. These are illustrated 
iu Fig. III. Higher order diagrams with similar 
qualifications constitute the rest ol I *. It is now clear 
that I through I" can include diagrams in which the 
projectile particle interacts with more Ilian oik* target 
particle through meson exchange. We therefore sc|Kiratc 
oul iIr* processes in which only one target |«ariicle is in- 
volved and w rite 


» •- 1 i 1*,., • •• • 2\ • 

# * • 

A.= V A„ I V A„ I • . '»«’> 

i i i >.i 

where K„, K m , etc. are the emit ribui ions mini the mtil- 
limcson exchange irreducible diagrams in wlucli mil) one 
target |Karlicle is involved. A are ll.c iiiiilliuicsoii ex- 
change irreducible diagrams with more than one target 
|iariiclc interacting with the projectile. W-. Iiavc Jcmon- 
st rap'd, lo the fourth order, that these terms liar., K., and 
K u I are very small compared to »■,. It was shown in Reis. 
12 and 13 that higher order terms such as K, lt arc also 
small when the scattering takes place at threshold. We 
believe that the higher order A”s will also !«• small above 
threshold. 

We note that it is possible to ileliiK separate I matrices 
with deriving terms K. r A’„,, e«e„ and derive a multiple 
scattering scries. In tit;*.: approach the resultant multiple 
scattering series will have a whole sv*t of » uiatric *s. Wc 
want to formulaic our multiple scattering scries iu terms 
of a new l matrix whose kernel is the sum of irrcducihlc 
diugrams involving only one target particle, i.v.. 
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|-|G. IU. Diagrams which contribute lo I' are shown to 
fourth order IIr- first term is the one meson evetunge term. 
The rest of the diagrams are the irreducible it tag rams as defined 
in the text. Ik dolled circle on a line indwales lb >1 the dia- 
gram is to he calculated without the on-shell eoninbution tor 
the projectile Hargett. ‘ITiese irreducible diagrams as a whole 
are defined as I'* in |q. tJ.4l. 


f ( I k", I * - • . As pointed out above. iIr* terms K m . etc. 
are small compared to i\ hut wc will include llnse terms 
in iIr* keriR-l of llie ir*w f matrix for the sake of complete- 
ness. Tlie terms given by A, involve more than one target 
particle and they will he treated in surli a manner that 
they would appear as higher order cttrrcclions in the 
resultant multiple scattering series. 

In lilt following we will show how* a multiple scattering 
series lor the optical |totenlial that ctMTespomls to lap 
13.31 can he ••biaincd. Ily employing the pntjeciHui 
operator nRtliod. we can rewrite l ; .q. 13.31 as coupled 
equations: 

f 0 » Cut;' j‘i u i fru* „ru . (3.7a> 

u r i i\k;; „qu r i i'Qo; „c ;>r . i3.7i»i 

Here (/ is our optical potential o|*cralor ami we seek a 
multiple scattering series e.xpressimi lor this operator. It 
should be noted that F’qs. 1 3-7* are ilircc-dimcnsional 
equations. llie first •me, l-.q. (3.7at. is the clleclive one- 
bod) equation I’m tin- projectile, and for a fermion projec- 
tile it becomes the fixed energy Dirac equation. 

Next, as in llie NR theory, we write V V, 1/ and ob- 
tain 

V U, £tf, 1 • • • A,» 

f » 

t Vtf. I 3 • 1,1 G*f „t', . 0.1 0 

'■i 

In the altovc equation, S' I »» the following, the projec- 
tion operators /* and Q art aippressed. It is to be under- 
stood that /‘goes with (i, 1 „ and Q goes with f«*’ Add- 
ing and siibiraeii.ig the qtiauliiy V, tr, I k„ < *** Igff, 
ami dropping the sum over f gives 

f/ if 1 , I A’„ I * * H if’, I I • • • Igf', 

+ Sir-, I A’., » ••• MG' K 

i 

• a, t v a,g; ,//, . 


t.i.-ii 
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covariant multiple scattering series for elastic . . . 
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where we have introduced a new propagator g whose 
properties are not specified at this stage Taking the 
second term on the right-hand side (RHSI of Eq. (3.91 to 
the left-hand side (LMS) and operating from the left with 
the inverse of [ I — (i», +R U + • ■ - |g] gives 

0,=7,+ |1+ 2g: o 0, | . 

where we have defined t, and r, a* 

T i =id i +K ii + ••• ) + (?,+£'* + • • fcF (3.MS 

and 

+(i»| + — )g9 f 

= A j +(t>, + + — IgA, + ••• 

^Ai+^gAi. (3.1 II 

Resumming over the index i gives 

0 - 25 + 2 ^: »-*>$+ 2 

• t • • # 

+ 2 p.|‘+2^»0>|j • »i2» 

Equation (3.121 is our multiple scattering series for tit- 
optical potential, and it is the second major result of this 
paper. It should be compared with the NR analogue Eq. 
(2.19). 

The first term of this series can he interpreted as the 
single scattering term for the optical potential and it is 
given by a r matrix driven by a kerne* which is liic sum of 
all irreducible diagrams involving only tnc target parti- 
cle. The second term is the propagator correction term 
which obviously depends on our choice oT the propagator 
g and its NR analog is the second term of liq. (2.19). The 
third term on the KHS of (3.12) is liic relativistic analog 
of the third term of Eq. (2.19'r and corresponds to multi- 
ple scattering corrections and they arc directly related to 
two, three, etc., particle correlations and can be assumed 
to he small in the first approximation. 

The last term includes iterations with I lie irreducible 
ding v.hich have more than one target particle in- 
vol' the definition of anil A, ) and it docs not 

con . in liic single scattering processes. These terms 
hav. „o NK analog and they originate from the nntl- 
ti meson exchange irreducible diagrams. Wc note that it 
is also possible to recombine this last term with the third 
term by separating the pieces according to the number of 
target particles involved. l*or example the first contribu- 
tion to the double .scattering would lx- £, ( , K ir It 
should he noted that because of the conspiracy among Ihe 
diagrams, these terms arc quite small and will not carry 
the same weight as the true double scattering term such 
as T,GZ in f r but they would constitute an in'-.-rcsling 
correction to the double scattering term. 

Keeping the first term only gives a single scattering ap- 
proximation for the optical potential. T he propagator# 
has not yet been specified. In principle, one could use 
any convenient i matrix for the i operator in Eq. (3.12) as 
long as wc are willing to incorporate the corrections 


represented by (he rest of the terms in Eq. (3.12). In 
practice one usually kci ys only ihe first term of the series 
and the judicious choice ofg is then essential. 

Under normal conditions, the second term gives the 
largest correction to the single scattering approximation, 
and we therefore should pay the greatest attention to this 
term. Wc would like to choose our propagator g so that 
this correct ion is minimal. This ran he accomplished by 
choosing the propagator e as shown in Fig. II. In this 
figure, both Ihe heavy A - I cluster and the projectile are 
kept on their mass sh i." In I lie medium energy range 
the terms riprrx -filed by the sum J , n Qt,- are dom- 
inated by the one nucleon knockout term and our choke 
of g described above would exactly cancel these dominant 
inclaslk contributions and ensure that they are exactly 
accounted for in Ihe t matrix itself given by Eq. 13.10). 
Restricting *he A — I cluster to the mass shell ensures 
cluster separability oT the remaining two-nuclcou sys- 
tcm. ,? 

With this choke of g, Eq. (3.10) for t, in the A7V sub- 
space reduces to the one panklc on mass shell (spectator) 
equation previously introduced by one oT us. 1 * A feature 
of lit is equation, discussed in Refs. X and 13. is that the 
multi meson irreducible diagrams lf„. K m . etc., are ex- 
pected to be small compared to and can be neglected. 
To be specific, after neglecting these terms the T t of Eq. 
(3.101 liecoincs 

/,=?, fi’.gt i ■ (3.10a) 

The projection o. /' onto this suhspacc will he denoted 
by I,. The only difference between /, and the free two- 
body l matrix is the shift in the total energy of the two- 
hixly suhspace due to tin- motion of the free A - I clus- 
ter. In analogy with the NK theory'- this choice of c can 
he viewed xs the “"impulse approximation" choice of g in 
our theory. The spectator Eq. (3. lift) is shown diagram- 
mat kally in Fig. 12. 

We conclude that the most appropriate I matrix ;» be 
used in the opt ice. I potential should he calculated from a 
covariant three-dimensional equation for two pankles in 
which one particle is kept on its mass shell. This choice 
will minimize ihe leading correction to liic multiple 
scattering series Eq. 13.12). 

The Iasi slep is lo carry out the nccissarv projections 



A 

I'll *. It. the uptime, it choice ol' ihe propagator v of li|. 
13.101. Ihe projti'lili- ami die A I cluster are both on the 
mass shell, indicated by a eross. I lit', choice of g mimmi/es tlx- 
leading eorreetiou term {the second term of Eq. I.V 1 2)|. 
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I Hi 12. The tact puiki Ml >Mli .xjujtMMi tv 

tlkv* lii iffj — i iIkiM) The aiiwo an a fine wtj* ilui iht 
(uniti- «v an ilk mass sk-H 


la obtain final equations Air elastic veal Icrme in ihc im- 
pulse affrouaiatKitt. Tit? dalk 'Olicnnf ampbink- n 
#’ - tlV. nhn.il viinfitv the iqua'iori 

r-y 4 vg' jr . 13.131 

«bcri V — PUP. ijjcaikta 1 3-131 r\ just Ike projveinwi <d 

laj. 1 3.7*1. In I hr smgk* sca«!cnng appro smtal mo. 

t/=Yf,. 13-141 

« 

where V - tl’t'. I tfualian 1 3.131 1 \ pra Ibr pf^vitm of 
E*M3 7 j|. i*< iht single vallitmt affriivniulMR. 

‘ *,*?, 13.1UW 

Using hath <3. Ivfcn and tJJUbl we obtain an alter nalivv 
iquatwm for 

t, ~r, ♦ r sf, l rgi,g?, . *3 I5« 

If wo define r,ll ( • lo tv iht* project khi i»f rtf, I .into the 
urinfnrc «4 states connected lo g. iht firsti»rder optical 
potential in ih, impulse ap-»ro'.tiHatwiw II A I «\ finally « At- 
tained a. 

VrfP 

- Vlff/ 1 ff.sr.r I Fr.Kt.sr. P* . <3 Ibal 

•hie 

t, - i* t 1 1 \gt 4 e t 1 ij>r. . i lib! 

Note ih.n lij. ( VltJ.i in cqwultm to th.- mo halj »t)iu- 
Inin with the project ik- on-shell. a\ divnU-il ahmt ami 
illuMratcd in lig. 12 ami i3. I$i tclk Imw all loer lip ol 
■hlN two-hody I matrix ate c»!rapolali-d «*ll nIi.I 1 lot MV n> 
the ojrtteal potiiili.il Note that I ico all tool kgs otl- 
nIvII land iiitlmkN a delta function hi the I 1 Njuvljlm 



I Mi IV Flos firm t hpcvii'i 1 Ik- <|tu.tiaieiv i-ijaitnii 
tilled I h* tutli ntT shell I mains » i> sji. -m « hi an njvn mat. 
I W itu WiS mat >s the speit ji.h t maim ol la) « t lib’ I lu- 
ll rsi icon on the K||S is tin- lulli isl slu-fl iiisuhi <<i one m.-vm 
tvluMf caulrihutmn us.,1 in I «| i t. tltai 


coordinates!. hut is identical to I, if one particle is un- 
shed tn the initial and final stale land the deha funetHm in 
the t I coordtuates ts drofvjvdt l-’nri htiHort. no Ittr- 
Iber eijualton must h’ solved to > At a in it: it is obtained 
directly front I, by quadrature. Kq I). Ihal. Equation 
1 Vital is illustrated in Pig. 13; its fourth-order cortribu- 
litw uras already enoumtefed in one of the terms in Ftg- 
v 


IV. IMStVSSION AN1) (UNCt.lINMJN 

The idea that the projectile-target scallenng amplitude 
is given by a relativistic equation where the kernel Of Ihe 
equal mo is the sum of certain diagrams was previously 
introduced in Ref. 17. However, these authors did not 
discuss nvcsoR exchange, nor dnl they derive any drtaik 
aboul the multiple seatiering series or the appropriate 
choice of relativistic equation and Ivro-hody amplitude. 

In this paper we have considered a euvartant formalism 
for pc< ijrr iik--tai get scattering in the eiMcil of meson ex- 
change. ami have shown that a multiple scattering scries 
for the optical ptatcniial can he denied. We d«» not claim 
that we have derived an RMSf starting from a field 
ihcoirlir'al i apiangian. hut we th* claim that we have de- 
rived a multiple scattering theory in a covariant manner, 
r.very step of inr derivation is manifestly covatiani and 
the end res. It. the t matrix associated with the impulse 
approximation optical potential. must ah** l*c calcubted 
from a relativist tcally covariant equal hut 

In cumTuston. we will restate what we have accom- 
plished mi the. pa|»rr. In the contest id - meson exchange 
we liave derived a invariant equal ton for the pritjectilc- 
nuek-us r ntairis lq l3 .31. I"hts equation uas thrived by 
cotisKkung the cancellalMins between the hos and 
crusscJ-lwi.x diagrams ami we Itave also shown t«* >» the 
spin mus singutai tlies can he avoided. We I hen derived a 
miiliiflc scattering series for the optical piHcnital and 
showed. in ihe impulse approsnnatiim. that ttie r mains 
associated with the optical potential is to he cakrulatcd 
from a relativistic thrcv-dinwnsMtiul equalttm in which 
one parin'!-.- is kept twi Hs mass shell. \Vc also described 
how llie fully off- slid I cxlrnstoo id" ilus i mains liq. 
l.l.Ufal eau be ci.Vululcd trow* a quailtaluie. I’q. t.Utu). 
We entphasi/e that our development leads t«» a precise 
definition of th- t minis to he used in the impulse ap- 
prosutiatioii «d the first -otilev opt teal potential. This is 
lie principal Jifli-icncc hetwrvw our result ami the RIA 
as ciMtmtonly used. 1 he r mains is to lie obtained from a 
one particle on mass shell equation. Hence intermediate 
stales with both nucksms in negative energy Mules eau.a*t 
■svttr. eseepl at the "end points." as tllusir .ietl in rig. 
1 1 lli.s result is obtainnl from a careful analysis of 
meson cvch.ingc diagrams, and seems to hi- tin- iim*si ap- 
pr.i|«rtatc lor tiie problem ol’ clastic nuehsMt-iiueknis 
-eatlcfing. Nttmcncal tests support this afq-roach. It lias 
hceit h*und that lh-‘ conlrihttltmts Itoiti rli tiinels m which 
kith iiiklisMis arc m* •attic energy -.tatev are iMghgi- 
hk‘. ’ I Ik - amp hi in hblivl ft om I ij fl ltd»> have* 

hivn used in au an., sis ol /•- '‘X .t cl.istw scallenng, 14 ' 
ami excel il agree at with espetimenial data has Inn 
obtained I >ifl- icn .-. between »-*«rs amt tlm «*!' I jon and 
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cofwijwcr mmmvE, scAftaiws setts for suulnc , . . 


Waito* 1 wwr« toaMe. hmt not targe. 

Smee **nr CfsHtn ter impnhse appro*»mal*t*n optical po- 
tential is dcrmst from a nwhiple icaitattog theory. it i* 
fmssMe to mate systematic eomxftMtns to the Irsl-oi-er 
taricohlicMs. We first mMmI to calotiate the four fcg off- 
todtl 1 mail** fro** the <qpttorat»re erpatto* O.l to) and 
then etahtoe other cwwikw. Por csamplc the tlouMt 


Foijtiy we .pitot opt that we hawe «n*t eoosideral the 
problem rf antisymmetry between the fwajpetife and the 
target particles, nor the self-snwtotent treatment of the 
4 -hWy tatpel stale. We haw mhm ifiwred the cwmpfer 
firms dne to spin. It is wry likely that the ppepskBe* 
target ardwjwiictiy cam he ehrsely approximated hy the 
Taiedto Watson* ptrscnprkw used in NX cahrniatinns. 
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Relativistic p- m Cn Mattie scattering observables arc .ilnlilnl using relativistic NN amplitude* 
obtained (rani the solution el a two-body relativistic cquiioe Results at 200 McV are p rora t ed 
lor three sets of NN laftaafa. two with pate pa eada v cc tot c cap lin g for the pion and one with a 
25% admixture of psewd oa cat ar c o apKag Differences h a w een the predictions of these three mod- 
els provide a atntarc of the anccrtaiatjr at const r ucting Dirac optical potentials from NN ampli- 
fades. 


It is wdl known that the relativistic impulse approxima- 
tion (RIA) gives s good descript ion of medians energy p- 
nackvs scattering observables. 1 In this approach the 
scattering of the proton is described by a Dirac equation 
with an optical potential (/of the form 

( 1 ) 

where the p superscripts (known as p spin) are positive 
(negative) for nucleons in a positive (negative) energy 
stale. This equation therefore expresses the optical poten- 
tial as a folding of the NN T matrix describing scattering 
of two initial nucleons in p-sptn states ip. pi) to two final 
nucleons in p-spin states ip'. pi) by a nucleon density dis- 
tribution p which depends on the p spin of the initial aad 
the final bound nucleon. 

Previous work 1 has clearly established that the success 
of the RIA depends strongly on the If 1 **"* and £/*"'♦* 
amplitudes, which in turn depend on knowledge of T- 
matrii elements in which at least one particle in the initial 
or f nal state is in the negative energy state. (In this paper 
these amplitudes will be referred to collectively as ’'nega- 
tive p-spin amplitudes.”) Furthermore, (he size of these 
negative p-spin amplitudes at low energy has been found 
to be very sensitive to the amount of y* coupling used for 
the pion. It has come to be accepted that pure y % y m cou- 
pling is required to give a satisfactory description of p- 
nuc<cus scattering at low energies. In this paper we report 
results which show that a surprisingly sizable admixture 
(25%) of y* coupling will still give good agreement with 


p-* 8 Ca data at 200 McV. provided the NN T matrix has 
been obtained from a relativistic meson exchange calcula- 
tion which fits the obs er ved NN data. This 25% admix- 
ture seems to us quite large; in the recent fit to the NN T 
matrix used in this paper aad repotted in Ref. 4, no larger 
value could be found, aad Flcisencr aad Tjou* found some 
litre ago that h was impossible to fit NN data with pure 
(100%) r s coopting. 

Before describing our calculation in detail, we briefly 
review previous work, focusing on results which depend on 
the form of the off -shell pion coupling. In the original cal- 
culations, 1 the values of the negative p-sptn amplitudes 
were inferred by expanding T in terms of five Fermi co- 
variants, and fitting these to the on-shell T l+M *’d*ta s 
It was found that this procedure is ambiguous and sensi- 
tive to the choice of the five covariants.* For example, us- 
ing the covariant 7^ [which corresponds to assuming the 
pion coupling b ps eu dos ca lar (PS)) or — (7*9-7) 1 
*(7*9- 7' )j/4m I [corresponding to the assumption that 
the pion coupling b pseudovector (PV)] gives identical re- 
sults in the (++,++) sector, but their extrapolations to 
the negative p-spin sectors are very different, resulting in 
generally poor fits at low energy (200 McV) when the 
7*7? covariant b used. One finds divergent energy depen- 
dence of scalar and vector potentials as energy decreases 
and consequently divergent virtual pair contributions to 
proton scattering when the yfy} covariant b used. To 
overcome thb problem Tjon and Wallace 7 adopted a 
meson exchange model of the NN interaction with pure 
PV ' pling for the pion as a basis for predicting all 
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matrix dements ia Eq. (1). Use of a dynamical 
■odd which describes the AW observables is ■ qualitative 
way over Ik 0-1000 MeV range was thoaght to be the 
best way to mtaimirc the ambiguity ia pcedictioa of the 
optical potcatiaL The optical pot e a ti a l based oa a com- 
plete *et of —fK»«t>« (referred to as IA2) pro d uces good 
a gr e emen t with e xperi m en tal resahs over a wide range cf 
projectile energi es from 200 to 100 MeV.* The divergent 
beh a vior of the potcaliaK at lower energies is abseat ia 
this case. 

la tins paper, we stady the sensitivities of the /Miadeas 
scattering observables to realistic variations ia the amoaat 
of y* vs fV coapliag of the pioa. We coatiaae to 
demand that oar NN T matrix be the solatioa of a relativ- 
istic eqaatioa with a meson exchaage kerad which has 
beca adjusted to give a qaaatitativdy accurate fit to the 
JVffdatx Two aew sets of relativistic NN amplitades are 
obtained from the solatioa of the relativistic eqaatioa ia 
which one particle is oa shelL* To easare that the resell- 
ing amplitades satisfy the Paali principle, the oac-bosoa 
exchange (OBE) kernels ased ia this ca leal a t io a were ex- 
plicitly aalisymactrixed. 4 - 1 * Hence the foar classes of 
ampUiac'es T**~**\ T l ~***\ T l ***~\ and 
T a *- ' ’ caa all be obtained through aatisymmetry. or 
time reversal invariance, from T l *~‘*** which we will 
refer to simply as 7* <_ * . and the amplitades T l * 

and can similarly be 

obtained from refereed to as The 

amplitades and ' are all taken to be 

zero. Finally, the kN coopting ased in these new solutions 
is a mixed coapiing of the form 

Xr i +( i_ X )jiazi. 

The parameter X varies the mix of pseadoscatar and psca- 
dovcctcr coapliag, and is defined so that the on -shell am- 
plitade is independent of X. When X is anity the coupling 
is purely pseadoscatar and when it is zero the coupling be- 
comes pare pscadovector. 

Two OBE models have been found which fit the NN 
data equally well, but which have significantly different 
amplitudes. 4 Ia model 1, only the four mesons a. 
or. m. and p are used. This is the minimal number needed 
to represent the long-, medium-, and short-range nuclear 
forces, and a very good fit to the positive energy NN am- 
plitudes is obtained when the parameter X has the value 
0.2S, which is 2Z% pseudoscalar end 75% pseudowedor 
(the fit tightly constraints the value of X). In another 
OBE model, model 2, the m N coupling is constrained to be 
pure pseudovector (X — 0) consistent with pair suppression 
and chiral symmetry. In order to fit the NN data equally 
well, two extra mesons. 6 and i), must be included. (The S 
meson is needed to get the correct splitting between '5o 
and J 5| central terms, which emerges automatically when 
X —0.25.) These two models allow us to explore the sensi- 
tivity to the amount of pj cu dos c alar coupling one may use 
and still obtain a good fit to the NN observables. They 
both differ significantly from the model used by Tjon and 
Wallace. 

The results for the polarized p-*C a elastic scattering at 
200 MeV obtained from the two models are shown in Kg. 


I. C al cu la t ions are baaed on the (A2 formalism of Ref. t. 
and we used the relativistic 4 ^Ca densities s u p p li e d by 
Horowitz and Sen*. 12 Su r pri sing ly, both models give a 
reasonable description for the p nuc l eus observables, and 
it caa be seen that the mixed coopting model gives superi- 
or results over the pure paendovector coupling case. How- 
ever. since the integral in (l) has only been evaluated in 
the tp approximation (in which the T matrix is evaluated 
oa shdl and factored out of the integral) and other effects 
sack as Pauli blocking and vacuum pola riz a ti on have not 
been included, it cannot be coa d adcd that the mixed cou- 
pling case will continue to give the best results after these 
effects are taken into accoun t . 

Wc would like to emphasize that the relativistic NN 
amplitude, used were the results of dyna mic al calcula- 
tions based on a relativistic equation, in which aB parame- 
ters were fixed by the NN data. No further adjustments 
ia any parameter were made ia calculating the p- m C a ob- 
servables. Calculations have also beca performed at 300 
and 500 MeV and the predictions agree with the data as 
wdU as for the case of 200 MeV. 

la order to isolate other possible mod el d e p en d en c e 
arising from negative-energy c omp one n ts of T. we com- 
pare models I and 2 above with a calculation of Tjon and 
Wallace 7 ia Fig. 2. In this case we standardize the com- 
parisons by replacing the amplitudes ia each 

case by the on -shell amplitudes determ in ed by Arndt rt 
a/. 11 First note that the two pseudovector models (model 
2 and Ijoa- Wallace) are ia dose a greemen t in spite of the 



0«m(deg) 


FIG. I. Predictions for ^-*Ci observables for the X— 0.25 
model I (solid line) and X —0 model 2 (dashed line) described in 
the text. The tun are dau from Ref. 1 1. 
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FIG. 2. Results for ike p - m Ca ohscrrabks with theoretical 
T*** *** amplitudes replaced Ly the Arndt amplitudes. aad 
theoretical amplitudes ia other p-spia sectors left unchanged. 
Model I (solid line), model 2 (dashed line), aad the calculation 
by Tjoa aad Wallace (dotted tine) ate shown. The data arc as 
ia Fig. I. 


fact that the Tjon-Wallacc model uses a Blankenbcdcr- 
Sogar reduction of the Bcthc-Salpcter equation which in- 
cludes ( ) channels (in which both of the initial or 

fint! particles have negative p spin) in the NN sector, and 
additional N& channels. 14 both of which were neglected in 
model 2. We have checked that the differences shown in 


Fig. 2 are nor due to the additional ( ) channels. 

Thus the sensitivity of Dirac results to the heavy mesoos 
aad the pr esen ce or absence of NA and AA ch a nn els or the 

( ) states of the NN channel seems to be minimal. 

Next, note that the larger differences between models 1 
and 2 attributable to the form of pioa coupling are sub- 
stantially unchanged, in this standardized comparison, 
even though the depth of the oscillation in A r has in- 
creased s om e w h at for both models. This shows that the p 
nucleus results are sensitive to fine details in the 
T l ++ +* } amplitude, but that the differences between 
rr.odch I and 2 cannot be attributed to these sensitivities. 

in summary, we have shown that the differences in the 
predictions of models I and 2. both of which lit the ou- 
shdl AW data very wdl, are due to the differences in their 
admixture of y 5 coupling, which cannot be uniquely deter- 
mined by the on -shell data. While this model dependence 
a significant it is surprisingly small c on s id eri ng the large 
admixture (25%) of y* coupling required by model I. 

Finally, it is amusing that the simpler model I. with the 
exchange of only four meso ns and a 25% admixture of y* 
coupling for the pioa. fits the ob ser v ab le s as well as it 
does. This result suggests either that some degree of pair 
nonsappressmn on the Born level may be allowed, or that 
the a counter terms required to control the y* part of the 
pion coupling may already be included as part of the phe- 
nomenological a exchange potential used in the NN mod- 
els. In view of the success here, it may be worth examin- 
ing the results of such mixed coupling models, both at 
lower energies where the sensitivity to mN coupling is 
larger, and in other reactions, such as electromagnetic 
p roc esses or process e s involving pion production or ab- 
sorption. 
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tioas ate important ia rtrirnaiaiag the shape of the theoretical aagaiar d utri ha li oa s Wc co at- 
pare the theoretical calcalalioas to a* data acd Bad that oar resaks are siarabr to the shape of 
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There continues to be interest ia finding a sensitivity to 
collective or coherent mechanisms in subthresbold heavy* 
ioa peon production. ,J This interest is enhanced by the 
experimental ability to measure neutral ptoos J ~ J which 
lends to data that are not complicated by final-state 
Coulomb effects on the outgoing pratts. The present mod- 
el interaction has been described briefly in previous 
works*' 1 where pioa spectra were calculated over a range 
of incident energies from subthreshold to relativistic 
values for **0 on ,2 C; however, those calculations used 
only the lowest -order (*— 0) angular momentum mul- 
tipole. We have now included all orders and disc o v er e d 
that the multipole expansion is highly convergent in the 
subthresbold region. For the first time with this model wc 
calculate the angular distributions for coherent *° produc- 
tion in C+C collisions at 85 MeV/nucleon at three 
different energy cuts. 

The Lorcnlz-in variant differential cross section calcu- 
lated in the projectile rest frame for a given target isos pin 
AflT) is given by 


d i o 
d'pJE . 


2s E.y' d 
hv (2 *M* dE F 


JdVr-lCnl 1 . 


(I) 


where o is the incident speed of the target as seen in the 
projectile rest frame. E m is the total pion energy, and V b 
the quantization volume which cancels an identical term 
coming' from the second -order amplitude C«. The total 
energy of the final stale E f is the sum of energies from the 
pion, recoil projectile, acd excited target of momentum 
pj~. Making the high momentum, forward scattering ap- 
proximation that p r . »/*„*• where p.y b the momentum 
of the recoil projectile plus pion as seen from the projectile 
rest frame, then d/dEF^ M (Erlp T 'C 1 )dfdp r .. Further- 
more. the scattering amplitude b assumed to be approxi- 
mately independent of the target solid angle. These as- 
sumptions are consistent with peripheral collisions of lew 
recoil and simplify the integrations over the target 


momentum pj*. The sums over all target bospins ti(T) 
of the differential cross section are then taken. The 
second-order amplitude 


Cn 


2XF| Kj | AfX/V | K, | /) 

AT 

E.+UMc'-H^'+irjl * 


(2) 


where | /) refers to the initial state of target, projectile, 
and the kinetic energy of the target. | N) b the intermedi- 
ate state where the target b excited to a •/***! *. T“1 
spin-, isospirt -coherent, bobar-analog, nodear state, and 
the projectile b excited to a similar state except that a nu- 
cleon has been excited to a A(3,3) isobar of mass 
Af«r z -I232 MeV and width r*-H5 McV. The final 
state | F) contains the pion that decays from the A aad it 
b assumed that the projectile returns to its ground state. 
Also, it b assumed that the decay of the A from the pro- 
jectile does not influence the target after the A b created 
so that the center -of -mass plane waves and the internal 
target wave functions b et we en the intermediate and final 
slates are orthogonal. A non relativistic Brcit-Wigner 
propagator for the A resonance b used where the energy 
dependence of the width* has been neglected since wc are 
concerned here mainly with angular distributions. 

The matrix dements for A formation and decay in Eq. 
(2) are evaluated using the techniques of second quantiza- 
tion and involves a great deal of Racah algebra. The for- 
mation matrix (N | | / > which excites a particle-hole 

stale in the target and creates a A-hole state in the projec- 
tile b reduced by Wick’s theorem to matrix elements of an 
interaction between particle hole and scattering states. In 
order to do an estimate calculation, a separable model A- 
formation interaction was chosen to be 

pi “2oo8(r)g((*, )g(£, )(S,o,)(T,- t f ) 

*ZC*(ft,)-C*(ft,), (3) 
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where g(r), {({,), mad f (& ) are Gaussian shapes which 
are functions of the coordinate magnitudes of the 
nucleus -nucleus center of mass, projectile nucleon, and 
target nucleon, respectively. This form was motivated by 
the work of Gaatde, Kammuri. and Ostcrfdd,* but 
simplified to take into account target Lorenlz contraction. 
The main features of this modd interaction are the spin 
and bospin flip where the target is cxdtcd to a giant reso- 
nance as well as the projectile with the added feature that 
a A is created in the projectile through the transition spin 
and bospta. The scalar product of spherical tensors of 
rank k allows for an exchange of orbital angular momen- 
tum #kk betwee n the projectile nucleon and target center 
of mass The target nucleons, although spin flipped, are 
left in their spatial ground state so that the matrix cle- 
ment in the target can be approximated by a Lorcntz- 
coatracted form factor. This model rep r esen ts a 
compromise where orbital rngular momentum changes 
are neglected in the target single-part idc states, but al- 
lows for a simple handling of the target Loreatz contrac- 
tion as described previously. 4 The strength of the interac- 
tion is givea by which is initially determined from two- 
body scattering and the factor 2 has been included since 


the strength of the A N* vertex is approximately twice as 
large as the NNtt vertex where quark theory gives 
gut. “ (6>/2/5 )/**». 14 This modd represents a step for- 
ward compared to our previous modd." The calculation 
of the A -decay matrix dement V | Vi\N) has also beea 
d e scribed pr e vio us ly. 1 1 

The nuclear, particle-hole states are spin-orbit coupled 
to produce the j values of the particle and bole and then 
coupled to produce the total angular momentum of the 
nucleus. The isospins of the particle-hole state arc also 
coupled to produce the total isospin of the nucleus. Then, 
linear combinations of these states with particle-hole 
codficieals are taken to produce the total nuclear 
state and each matrix dement contains 9 / symbols after 
lengthy calculations in angular coupling are done. The 
panicle and hole states are generated from the three- 
dimensional harmonic oscillator and arc described in Ref. 
6 except that we now apply the modd to ,J C on IZ C. The 
sum over intermediate states i N) in Eq. (2) contains 
sums over panicle-hole states for both nuclei as well as a 
sum over angular momentum multipolcs k as given by Eq. 
(1) in Ref. 6 and described therein. This sum contains a 
product of terms which is gi v en by 


Tf~jt (Kr)g(r)r 2 drf~RlyOj t <.kJ)R'± (Ofrt [ q 0 q 


2 

Pk (cos© jr ) . 


(4) 


where the first integral is in the Born approximation with 
r espect to the target momentum transfer hK. and the 
second integral comes from the projectile A-holc decay 
matrix dement with respect to the outgoing pioo momen- 
tum hkf The p rop ert ies of the 3 j symbol 

U k /a] 

0 0 Oj 

restricts k to even/odd values if /*+/* —even/odd values 
and severely truncates the k sum through the triangular 
condition |/ 4 — /*| <l</*+li. The pion angular dis- 
tribution b contained in the Legendre polynomial of order 
k with respect to the correlation angle 0,x which is the 
angle between the pion and target momentum transfer. In 
previous work. 4 we only hand calculated the differential 
cron section for the k — 0 term and only calculated the 
square of Eq. (4) for k —0 plus Jfc — 2 to determine the rel- 
ative effect of the higher-order multipole. Presently, wc 
have a complete code which includes all k multipolcs and 
have calculated the differential cross section as a function 
of the correlation angle. After including higher A-hole 
stales, wc discovered that the main contributions to the 
cross section come from the Ipi/i.j/ivrA states and the 
lpj/j-holc state in U C (a) for the k — 0 and 2 multipole 
values. For example, when we included the A stales for 
the ld-2s shell which corresponds to k » 1,3 and * — I, re- 
spectively with the Ipyi-holc state, we found these contri- 
butions to be negligible. This occurs because the overlap 
between the functions in the first and second integral in 
Eq. (4) diminishes considerably and the value of the 3 j 


I 

symbol 

[/a * U 

[o 0 0 

drops. It appears that the k sum in Eq. (4) b highly con- 
vergent. Also, since wc assumed equal weight particlc- 
hole coefficients A^CP), they become smaller when more 
shell-model states are included. If these coefficients were 
calculated in the schematic model. 1 * wc would find that 
they arc inversely proportional to the particle-hole energy 
difference which cuts down the contributions from higher 
shell-model stales and justifies our truncation to lowest 
stales. We also included the 1/ shell and found that these 
results were also negligible. Therefore, the main contribu- 
tor to the nonisotropic shape of the differential cross sec- 
tions comes from Pi although the actual shape is deter- 
mined by the relative weighting between Po and Pj in the 
multipole sum. 

For the first time, we compared the results of this model 
to x° angular distributions in U C+ '*C collisions at 84 
McV/nuclcon transformed to the nudeon-nudeon, cen- 
ter-of-mass system.' 1 Since wc are considering small 
momentum transfers in the forward direction, wc assume 
that the correlation angle 0,* « approximately the pion 
angle measured from the forward direction 0,. For fixed 
values of the pion kinetic energy and angle (77^ .0.* 1 } 
in the nudeon-nudeon center of mass, we calculate the re- 
lativist -.cal !y transformed set (Tf,0,f) for pions produced 
by the projectile in its own frame to obtain the Lorentz- 
invariant differential cross section which is further 
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transformed to the noninva riant cross section in the 
center-of-mass system. A similar procedure is applied to 
oootribution coming from the target These two cross sec- 
tions are then added so that oar results contain the in- 
coherent addition of pions coming from projectile and tar- 
get; however, the A's are produced coherently in each nu- 
cleus. The results are shown in Figs. 1-3 for three 
different energy cuts. Since we are concerned primarily 
with shape fits, the calculations have been normalized to 
the data at cosB. “0.025 since that data point has an error 
bar shown in Fig. 3. We do not expect to fit the data in 
absolute value since the calculation is exclusive whereas 
the data is inclusive. It is interesting to note that the 
theoretical shape improves at the highest -energ y bin from 
100-150 MeV where the peak-to- valley ratio moves to- 
wards the data. Hus trend is compatible with the con- 
clusions of Braun-Munzinger and Stachd 1 where they 
point out that the low- en e rg y pions are probably emitted 
from a local hot spot, whereas the pions produced approxi- 
mately above 100 MeV cannot be fully understood by 
thermal models aiooc and may suggest the pre sen c e of a 
coherent mechanism in the production process. It is also 
interesting that our normalization factor drops by a factor 
of 10 from the lowest -en er gy cut to the highest perhaps in- 
dicating a conv er g en ce at the higher pion energies. It is 
also intriguing that for the lowest -energy cut. the results 
from the bremsstrablung model compared to similar 
data 14 shows a convex trend as do our results under the 
condition of incoherent ad "lion of projectile and target 
pions. whereas the bremsstrahlung calculations and our 
results improve at the higher-energy cuts. 

In conclusion, for the first lime, we include pion contri- 
butions coming from the target and projectile and find 
that these contributions are important in determining the 
theoret ic al shape of the r° angular distributions. Wc 
compare the present model with experimental a° angular 
distributions but do not expect to obtain absolute-value 
fits to the data because our calculation is exclusive 
whereas the data to which we compare arc inclusive; how- 
ever, we wanted to find out if the qualitative shape and cn- 





FIG. I. Theoretical x* production compared to (he data (Ref. 
12) of C+C collisions at 84 McV/nudcon for the 0 to SO MeV 
cal. 
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FIG. 2. Same as Fig. I except for the SO to 100 MeV cut. 


ergy trend of the theory bears some resemblance to exist- 
ing experimental angular distributions before embarking 
on the more complicated phase of this calculation which is 
to include the tensor term. We find that the theoretical 
shape displays a forward-backward peaking similar to the 
shape of the data at the higher pion energi es. The for- 
ward peaking in our model is due to pions coming predom- 
inantly from the projectile whereas the target predom- 
inantly contributes to the backward peaking as seen in the 
nuckon-nndcon center-of-mass system. Wc expect this 
qualitative result would also obtain when the tensor in- 
teraction is included. The calculations is also very sensi- 
tive to the shell-model information since the multipolc 
sum is very conv er ge n t and not much of a washing out of 
the shell-model signature occurs. This was not obvious 
from previous work 4 * 7 where hard calculations were done 
and it was not known if a washout of the angular distribu- 
tion would occur with the inclusion of higher-order mul- 
tipolc values. We again expect a similar result with in- 
clusion of the tensor term and now have a tentative ex- 
planation for the forward-backward peaking seen in the 
data. 
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FIG. 3. Same u Fig. I except for the 100 lo 150 MeV cut. 
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The model isobar-formation interaction, even though 
incomplete, leads to a complicated, microscopic, quan- 
tum-mechanical. many-body, angular momentum formal- 
ism. Our approach was to start with the simpler spin-spin 
term but include orbital-angular momentum exchanges 
be t we en projectile and target nucleons as well as the 
scattering of the relative nucleus-nucleus system and cal- 
culate angular distributions to discover the nature of the 
signature for coherent production. From Ref. S, in an ex- 
amination of the momentum transfer dependence q of the 
central and tensor interactions, it is seen that for low 
momentum transfers of 0.5 fni -1 . the central term 
dominates over the tensor term. In fact, at q “0. the ten- 
sor term is aero. At increasing values of q, the central 
term drops and the tensor term rises until at the critical 
value of q ** 3 fm ~ *, the tensor term dominates the cross 
section because the central term goes to zero. However, 


for peripheral collisions at sabthrcshold energies, it is like- 
ly that use lowest values of momentum transfer will be 
favored and that the central term might dominate these 
reactions. With this work, we are encouraged to include 
the more complicated tensor interaction and compare the 
roles of the central to the tensor term. In Ref. 2, where 
the tensor term was included, only total cross sections 
have been calculated for the subthreshold process using 
the impulse Fcshbach-Zabek 11 approach. In our work, we 
are attempting a more fundamental, microscopic, quan- 
tum-mechanical approach. We are presently developing a 
hcavy-ion calculation with a more realistic interaction 
model that includes the tensor term. 
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RcUtivnlic p-*Ca elastic scattering observables are calculated using four sets of relativistic NN 
amplitudes obtained from diSerent one-boson-eichangc (OBEI models. The first two sets are based 
upon a relativistic equation in which one panicle is on mass shell and the other two sets are ob- 
tained from a quasi potential reduction of the llrihe-Salpcter equation. Results at 200, 300. and 300 
MeV are presented for these amplitudes. Differences between the predictions of these models pro- 
vide a study of the uncertainty in constructing Dirac optical potentials from OUE-based NN ampli- 
tudes. 


I. INTRODUCTION 

In this work relativistic />- nucleus calculations based on 
four different sets of NN amplitudes obtained front rela- 
tivistic meson-exchange models are presented. The com- 
parison of the results based upon these different relativis- 
tic NN amplitudes is motivated by the fact that the main 
ingredient in the success of the Dirac approach to p- 
nucleus scattering is the virtual NN effects 1 which are im- 
plicit in the solution of the Dirac cqualion. For the posi- 
tive cnetgy component of the wave function wc can write 

lE~E f — U 1 * =0. 
where 

u^, = u' ie+E'-u r'ir ' . 

To a good approximation U ! ' is the nonrelativist ic 
potential and it is essentially determined by the NN on- 
shcll data and the density prolile of the target nucleus. 
When one uses meson-exchange models >o calculate the 
relativistic NN amplitudes the pure positive energy am- 
plitudes are fitted to the on-shell data and the amplitudes 
which belong to the pure relativistic sectors are thee 
determined dynamically. In principle all meson- 
exchange model calculations should provide the same 
pc.itivc energy ampli’udcs (since they arc lilted to the 
on-shell data) but the amplitudes in the other sectors 
would be different depending on the dynamical model 
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used. Therefore the ambiguity in the Dirac optical po- 
tential lies in 'he U ^ term, which is determined by am- 
plitudes other than the on-shell one. Since different 
dynamical models will give different amplitudes, except 
in the pure positive energy sector, it is important to study 
these differences and possible effects on the p-nucleus 
scattering observables which may arise because of these 
differences. 

In previous work,* preliminary analyses of relativistic 
p- 4 °C a elastic scattering at 200 MeV, using relativistic 
nucleon-nucleon (NN) amplitudes obtained from the 
solutions of two-body relativistic equations, were report- 
ed. In the present work the calculational methods are de- 
scribed and an extensive array of results for 200, 300. and 
50C MeV protons incident upon a '* u Ca target are present- 
ed. Calculations arc preformed with two sets of NN am- 
plitudes obtained from a relativistic wave equation in 
which one particle is on positive energy muss shell 1 and 
also with third and fourth sets of NN amplitudes based 
upon a qui'sipotcniial method 4-5 for solving the Bcihe- 
Salpeier equation. The use of the relativistic NN ampli- 
tudes obtained from a one particle on-mass shell cqualion 
is motivated by a recently formulated covariant multiple 
scattering theory. 6 In this formalism the authors argued 
that the NN amplitudes associated with the rclafivistic 
itnpulse approximation optical potential should be the 
one obtained from a one particle on-mass shell cqualion. 

The outline of the paper is as follows. In Sec. II wc 
present an overview of the optical potential methods used 
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lo compute the elastic scattering observables. A brief dis- 
cussion of the Tjon and Wallace method 7 for obtaining 
the invariant AW amplitudes is also included for com- 
pleteness and further discussions. In Sec. Ill the meson- 
exchange models used are described together with the as- 
sociated relativistic equations and the one-boson- 
exchange (OBE) parameters are tabulated. !n Sec. IV re- 
sults for elastic scattering observables are presented for 
all four sets of NN amplitudes. To test these results for 
possible model dependencies arising from the differing 
negative energy components of the r matrix, the on-shell 
amplitudes are replaced by the ones determined by Arndt 
er o/. ( and the resultant scattering observables computed 
and discussed. Finally, Sec. V summarizes our major re- 
sults and conclusions. 

II. THE DIRAC OPTICAL POTENTIAL 

It is well known that the relativistic impulse approxi- 
mation (R1A) is a useful method for predicting intermedi- 
ate energy />-nucleus scattering observables.’ In this ap- 
proach one begins with the fixed energy Dirac equation 

[yp~ tn ~ Olr)]4ir)=0 , (I) 

where the energy E is £=(p ! +m J l l/! , p is the on-shell 
momentum, Y-p—Y^.p"' — (£,ph and the repeated in- 
dices are summed over. The optical potential 0, obtained 
from the RIA, is written in momentum space as lu 

f)(p',p)= — ~Tr, 

X J JdkM(p,k -iq-*p’,k + {q)plk,q) J , 

( 2 ) 

where M is the Feynman NN amplitude and ji is the rela- 
tivistic nuclear density matrix for the target. In Eq. (2) 
particle 2 is the target nucleon and the trace is taken over 
its Dirac indices. Using the factorization approximation 
Eq. (2) yields 

0{p',p)=-{Tr 2 [M(p,-\q-^p‘Aq)p(q)) . (3) 

Therefore to construct 0 we require M and />. For this 
work involving *’Ca target, the relativistic densities of 
Horowitz and Serot" are used for />. The remaining 
problem is then to specify M. 

It is customary to separate the full NN amplitudes into 
a sum of 16 terms each representing a particular p-spin 
sector. They arc labeled by indices whe- each 

p is either + or — , and primed (unprimed) quantities 
denote final (initial) nucleon states. Tjon and Wallace 
have determined M for all 16 p-spin sectors by solving the 
Bethe-Salpeter equation using a quasipotential method 7,5 
with pure pseudovector coupling for the pion. In that 
particular formalism, the invariant NN amplitude is ex- 
panded as 

/>=--* = y . (4) 

x , , 
p 

where each p is summed over T and — , and k is a kine- 
matic factor [see Eq. (2.10) of Ref. 8]. In Eq. (4) the 


arc am p|it U( |es restricted to act only within one 
p-spin sector (hence there are 16 different sets of them) 
and the A*s are covariant projection operators given by 


A'(p>= 


2m 


(5) 


where p is + or — . Details for determining the various 

p-spin sector amplitudes are found in Refs. 7 and 

10. This complete calculational program is referred to as 
the 1A2 formalism. 10,12 


HI. RELATIVISTIC MESON-EXCHANCE MODELS 

In this section we briefly describe the relativistic 
meson-exchange models used in this work. In operator 
form the relativistic two-body equation can genetically be 
written as 

M — V+ VGM , (6) 

where M is the invariant amplitude for the two nucleons 
scattering process, V is the interaction determined by the 
OBE modci, and G is the relativistic two-body proftaga- 
tor. M satisfies the Pauli principle provided that V is 
suitably antisymmetrized. 

The first two OBE model amplitudes 17 used in this 
work are calculated from a one particle on-»nass shell 
equation and in this case the propagator G in Eq. (6) is 
the propagator with one particle on the mass shell. Fig- 
ure I depicts Eq. (6) for this case diagrammatically. The 
open boxes represent the combined OBE terms, properly 
antisymmetrized, for the mesons used. The horizontal 
lines represent the nucleons and a cross on a nucleon line 
indicates that it is on the positive energy mass shell. Ex- 
plicit antisymmetrization of the OBE kernel insures that 
the M matrix is correctly antisymmetrized under particle 
interchange. The first OBE model uses four mesons 
( *r, a,p, and co). This is the minimal number needed to 
represent the long, medium, and short range nuclear 
forces. The irN coupling used is a mixed coupling of the 
form 



FIG. I. Diagrammatical representation of one particle on- 
mass-shcll equation. A cross on a nucleon line indicates that 
the particle is on the positive energy mass shell. Open boxes 
stand for rhe anlisymmetrized kernel, and vertical dashed lines 
stand for combined v, a, p, and to exchange for model I and 
it, a, p, to, ?/, and a , exchange for model 2. 
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FIG. 2. Diagrammati il representation of Bethe-Satpctcr equation for coupled channel NS + AfAtAA scattering. The A stales 
are indicated by double horizontal lines. The dashed line indicates jr +p exchange: the wiggly line stands for combined «r. rj. r, b. p, 
and u exchange. Tlic single SN channel case can be obtained by leaving out the diagrams with A lines. 


Ar’+n-A)-^. 

2m 

where the parameter A. varies Ihe mix of the pscttdoscalar 
and pscudovcclor coupling, and is defined so that the on- 
shcll amplitude is independent of A. When A is unity the 
coupling is purely pseudosealar and when it is zero the 
coupling becomes pure pscudovcclor. Use of this cou- 
pling permits us to study the sensitivities of the scattering 
observables to relative variations in the amount of y f vs 
coupling of the pion. The best value of the mixing 
parameter A to fit the data was found to be 0.22 which 
corresponds to a 22% y 5 coupling admixture. 


In the second model, the irN coupling is constrained to 
be pure pseudovector ( A — 0 ) consistent with pair 
suppression and chiral symmetry. In order to fit the NS 
data equally well with this constraint, two additional 
mesons »/ and a , were added. The it, meson, a phenome- 
nological spin zero isospin 1 meson, is needed to obtain 
the correct splitting between the '5,, and '.V, central po- 
tentials wiien A =0. Note that neither of these OBE mod- 
els contain A/ A or AA channels but both models are able 
to fit the NN phase shifts very well up to 300 MeV. 

btnee in these models (I and 2) we require that one par- 
ticle is on the positive energy mass shell p-spin ampli- 
tudes labeled by ( ,p,p> ) ..nd I arc zero by 


TABLE I. 

Meson masses (MeV), cutoff masses 

(MeV). and 

TABLE II. 

Meson masses (MeV). culolT masse. 

(MeV), -*nd 

coupling constants for OBE models I and 2. 


coupling constants for OBE models 3 and 4. 



Ml Kill 1 

Model 2 


Model 3 

4 

« r j?;/4n- 

13.544 03 

13.357 57 

, 

~ ~ ----- ■ 

- *"■ 

m _ 

138.0 

138.0 

n g;/4ir 

14.2 

14.2 

3, 

0.225 57 

0.0 

m _ 

138.7 

138.7 

a gi /4b- 

5.513 22 

5.04720 

b K.i/4n- 

0.33 

0.33 


516 0 

514.0 


960.0 

960 0 

P ?r /4ir 

0.38291 

0.58686 

9 K',/4i r 

3.09 

3.09 

K r 

7.525 25 

6.149 20 


548.0 

548.0 

K 

1.0 

0.752 18 

P Kr /4rr 

0.43 

0.43 


760.0 

760.0 

K. 

763.0 

763.0 

to J?i/4ir 

9.85 106 

9.83054 

P Kl. /4ir 

6.8 

6.8 


0.142 59 

0.150 50 

< 

763.0 

763.0 


782.8 

782.8 

•o K'. . /4ir 

1 to 

12.0 

Oi 


0.325 93 

Hi M 

783.0 

783.0 




f c, /4ir 

7.65 

4.7 



573.0 




i 



m, 

570.0 

570.0 

V ir,/4tr 


6.40798 

/ii/4ir 

0.0 

0.35 

m. 


548.8 

fU /An 

0.0 

1.0 

Ay 

1610.0 

1685.0 

Aw 

1150.0 

11500 

^(Hexnn 

2135.0 

1830.0 

_Aw._ 

1070.0 

1070!) 
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definition. In addition, the ( 4 — . + 4 ). ( - +. 4- 4- ). 
( + +, + -). and( + +.- + ) amplitudes, which we will 
label simply as (—. + ), are obtained from the 
( + —, + + ) amplitudes using anti s y mm etry or time re- 
versal properties. Similarly. the amplitudes 
( + -. + -), ( + -.- + ).(-+. + -). and (-+.- + ). 
labeled as (—, — ), are obtained from ( + —, + — ) ampli- 
tudes. 

The third and fourth OBE model amplitudes*-’ that we 
employ in this work are calculated from a quasi potential 
equation based upon the Bethe-Salpetcr equation. In 'his 
case the propagator G of Eq. (6) is the quasi potential 
Bhuikcnbecler-Sogar propagator.*-’ The third model in- 
cludes only the NN channels.* In our calculations, the g ( 
coupling co n stant of model 3 has bee n increased slightly 
as compared to Ref. 4, in order to improve on the fits of 
the S-wave phases. It is needed because of the inclusion 
of the negative energy spinor state channels, which has a 
repulsive nature in the S-wave channels. In addition to 
the NN channel the fourth model also includes the N A 
and AA channels ' The scattering equations for these 
cases are shown diagrammatically in Fig. 2. The equa- 
tion for the case where there are no N A and AA channels 
can be obtained by leaving out the diagrams with the A 
lines. 

The meson masses, coupling constants, and cutoff 
jaws for these models are given in Tables I and I 1 . To 
see the goodness of fit to the on-shell NN data, the Wol- 
fenstein A and C amplitudes at 200 McV are shown in 
Fig. 3 for the first three OBE models together with the 
empirical Arndt amplitudes. It is to be noted that while 
the firs* two models resulted in reasonably close fits to 
each other, the third model differs quite visibly from the 
first two and also from the Arndt amplitudes. This gord- 
ness of fit to the on-shell data is also one of the possible 
ambiguities that can arise in the optical potential when 
the NN amplitudes obtained from OBE models are used. 
At 300 and 500 MeV, models 1. 2, and 4 give reasonable 
fits to the Arndt A and C amplitudes and they are shown 
in Fig*. 4 and 5. The fits of models l and 2 are generally 
less good than at 200 MeV. but it turns out that the fils to 
the NN amplitudes at 300 and 500 McV arc not as crucial 
to the success of the p- nucleus prediction 

IV. RESULTS 
Elastic scattering abscrvakles. 

The results for p-*“Ca elastic scattering observables arc 
presented in this section. Calculations arc based on the 
IA2 formalism of Refs. ^ and 10 and use the relativistic 
nuclear densities supplied by Horowitz and Scrol." 

figures 6-1 display results for elastic scattering ob- 
servables, obtained with the OBE models discused above, 
for proton encigics 200. 300, and 500 McV. In all three 
figures, results from the mixed coupling model 1 ' arc 
shown by a solid line, results from the pure pseudo vec I or 
model 11 are shown by a dashed line, and one of the two 
models of Ref. 4 is shown by a dotted line. In Figs. 6-8 
p-nuckrus calculations are done by using the theoretical 
amplitudes. For comparison purposes experimental data 


from Refs. 14. IS, and 16 are included for 200. 300, and 
500 MeV, respect i vel y . At 200 MeV, JVA and AA chan- 
nels are omitted in all three models. At this energy, mod- 
el 4, which is designed to study the effects of coupling to 
the A states at higher projectile energies, is not used in 
the calculation. All models give a reasonable description 
of the scattering observables, although model 3 (Ref. 4: 
with no NA or AA channels) s ee ms to describe the A r 
better in the forward angle region. Note that, although 
model 3 seems to describe the 200 MeV p-nuclcus data, it 
docs not reproduce the Wolfcnstein A and C very well as 
can he seen from Fig. 3. At this stage our conclusion is 
somewhat uncertain. It would he an interesting study in 
the future to compare the p-nucleus results of these mod- 
els with the same quality of fit to the NN on-shell ampli- 
tudes. These results also show that the two OBE models 
of Ref. 13 U=0 and A. =0.22) give similar results al- 
though model I (A. =0.22) seems to describe the experi- 
mental data somewhat better. 17 

For 300 and 500 MeV. the dotted lines are the results 
of model 4. which includes N A and AA channels, while 
the other two models (I and 2) remain unchanged. At 



FIG. 6. /v*'Ca scattering observables at 200 McV. Model I 
results are shown by a solid line, model 2 results by a dashed 
line, and model 3 results by a doited line. Note that in all three 
calculations NA and A A channels are omitted. Data arc from 
Ref. 14. 
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ihne energies, ail three models give reasonably good fils 
to the on-shell Arndt NN amplitudes, although there are 
some differences. From the 300 and 300 Me V results, it 
seem* that the experimental data do not appear to favor 
any one of the NN amplitudes over the other, although 
the polarization observables at 300 MeV seem to be 
somewhat better fit by models I and 2 in the region of the 
first minimum. We also notice that at these energies the 
differences between the A. =0.22 (model II and A. =0 
(model 21 appear to be very small suggesting that the p- 
nucleus scattering observables are not very sensitive to 
the relative admixture of y' coupling at these energies. 

As mentioned above, the best way to compare the p- 
nucleus results of these different models would be to have 
the saute goodness of fit to the NN on-shell data. Since it 
is impractical to do this at the current stage of develop- 
ment we proceed as follows. In order to separate the 
effects of the negative energy components from 
differences which may arise from the goodness of fit. we 
standardize the comparisons by replacing the 1 H . H I 
amplitudes in each model with the on-shell amplitudes 
determined hy Am' l i rf «/.‘ The results are displayed in 
Kip' ‘• — II. First note* that all three models are in reason- 
ably close agreement with the experimental data although 
the previous agreement with data and model 3 at 2111 



°em 


I IC. 7. scattering observables at .Ml MeV. Models I 

and 2 results are shown by solid and dashed lines, respectively. 
Model 4 Iwhreh contains \i and A A channels) results are 
show n by a dolled line Data are from Ref 15. 


MeV seems to be destroyed somewhat. This shows that 
at this energy '200 MeVl ^-nucleus results are very sensi- 
tive to the positive <wrjy sector of the NN amplitudes. 
Standardized calculations based on model 4 al 200 MeV 
iAfA and AA states included) have been shown in Ref. 12 
to give a very good description of the data. 

At 300 and 300 MeV, the mixed coupling model (mod- 
el II and model 4 are in very dose agreement in spile of 
the feet that model 4 uses a BUnkenbccler-Sugar reduc- 
tion of the licthc-Salpctcr equation which includes 
< — . — I channels (in which both of the initial and final 
panicles have negative p spin), and additional N A chan- 
nels, both of which are absent in models I and 2. There- 
fore the sensitivity of the results to the presence or ab- 
sence of/VA and AA channels or ( — . — I stales appears to 
he minimal al best, indicating that the /»- nucleus theory al 
high energy tat least up to 300 MeVl is not critically sen- 
sitive to the goodness of the til to the NN system. Com- 
paring Figs. 6-8 with Figs. 1-11. we also note that the 
difference between models I and 2. attributable to the 
form of the piun coupling, are substantially unchanged in 
the standardized comparison. 

For models I and 2 al 200 MeV. use of the Arndt am- 
plitudes appears to give a slight improvement in the 
agreement with the experimental data. In addition the 



FIG. X. p-*Va scattering observables al 5tt) MeV. The 
meaning ol I Ik- curves is iIk- same as m Fig. 5. I Tula are fr«*m 
Kef. id. 
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minima in A f appear deeper and the minima in Q appear 
shallower with the use of the Arndt amplitudes thereby 
demonstrating that the p-nudeus scattering obs e rv a bles 
are sensitive to the details of the (+ +.+ +) amplitudes, 
at this energy, even though differences be tw een model . 
and model 2 are not attributable to these sensitivities. 

V. SUMMARY AND DISCUSSION 

In this work we have presented p-nucleus scattering 
calculations based on four different sets of relativistic NN 
amplitudes. The first two sets of amplitudes are obtained 
from a three dimensional relativistic two-body equation 
in which one particle b kept on the positive energy mass 
shell- The difference betwee n these two models is that 
one uses a pure pseudovector wN coupling and the other 
uses a mixed coupling (22% pseudoscalar). The results 
for *00. 300, and 300 McV proton energies are presented. 
For all three energies both models include only the NN 
channels. Comparisons of the results of these two models 
show that, although both models give reasonable agree- 
ment with the experimental data, there are some definite 
model dependencies. This model dependency seems to 
lessen at higher proton energies. Since these two models 



0cm ( d eg) 


FIG. 9. p-®Ca scattering observables at 200 McV. The 
meaning of the curves is the same as in Fig. 4. but the positive 
energy amplitudes in all three models arc standardized to the 
Arndt amplitudes. 


reproduce the on-shell NN data equally wed, as shown in 
Fig. 3. we cannot attribute these differences in p-nucleus 
observ a bles In the ambiguity in fitting the on -shell data. 
This statement is also supported by the fact that the 
differences in p-audeus results seem relatively unchanged 
when the Arndt amplitudes are used tc standardize com- 
parisons. 

The third and fourth set of NN amplitudes are based 
ou the quasi potential reduction of the Bet be-Sal peter 
equation. The third model includes only the NN chan- 
neb while the fourth one includes both the NA and AA 
channels. Both models use pure pseudervector wN cou- 
pling. The p-auckus calculation at 200 McV was done 
with amplitudes obtained from model 3 (see Ref. 12 for 
one based on model 4) and at 300 and 300 McV model 4 
amplitudes w ere used. Since model 3 does not reproduce 
the NN oo -shell data at 200 MeV very well, there is some 
ambiguity in comparing its results to models I and 2 at 
this energy. When Arndt amplitudes are used lor stand- 
ardized comparison, the results of models I and 2 change 
in a minimal way while the model 3 results change quite 
visibly. This observation shows that p-nodcus observ- 
ables are very sensitive to the positive energy NN ampli- 
tudes at this energy. This in turn suggests that in order 



0cm(deg) 


FIG. 10. jhTj scattering observables at 300 MeV. The 
meaning of the curves is the same as in Fig. 5, but the positive 
energy amplitudes in all three models arc standardized to the 
Arndt amplitudes. 
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FIG. 11. p-*"Ca scattering observable* al 500 McV. The 
meaning of I be curve* K (he same a* in l-lg. 6. hul the positive 
energy amplitudes in all three models ate standardized to the 
Arndt amplitudes. 

to use relativistic A'.V ampliludes obtained from OBE 
models without ambiguity, it is essential to have a good 
quality hi to the on-shell data. Although this statement 
is quite obvious and seems unnecessary, it is important to 
note that p-nucleus scattering is one of the areas in which 
the pure relativistic amplitudes can be tested and these 
ampliludes are determined only dynamically (through 
lilting the on-shell ampliludes) once an OHE model is 
chosen. 

At 300 and 500 McV, model 4, which includes the del- 
ta isobar channels, and models I and 2 are compared. 
They all describe the experimental data very well. In par- 
ticular, the results for models I and 2 arc very close and 


there are no obvious model dependencies. Ftcepl for the 
fact that model 4 seems to give a very t'ecp first minima 
for A r and Q at 500 MeV, all th/cc models give very 
similar predictions for the p-nucleus scattering observ- 
ables. When the standardized comparison is made, the 
differences in the models disappear almost totally. These 
results suggest that at high energies (up to 500 MeV). p- 
nucletts theory is insensitive to the quality of fit to the NN 
ampliludes. including the presence or absence of N A and 
AA channels, and the ( — . — ) channels in the NN sector. 
Most of the effects of the A's in elastic p-nucleus scatter- 
ing may be described as giving rise to an attraction in the 
NN chaiittek at medium range and which are effectively 
simulated by a stronger epsilon coupling constant in the 
other imidd. 

In summary, we conclude that OBE models provide an 
unambiguous way of determining the relativistic NN am- 
pliludes provided that a good quality fit to the on-shell 
data is achieved. There is a definite model dependency in 
^nucleus results depending on the choice of wN coupling 
but this dependency seems to disappear as the incident 
proton energy becomes higher. It is well known that the 
original RIA. which uses the five Fermi covariants, gave 
unreasonably large scalar and vector components at low 
energies. This has been attributed to the large pair con- 
tributions which result from pseudoscalar wN coupling. 
Since one of the OBE models (model I) studied in this 
work contains 22% pseudoscalar coupling, it seems to 
suggest that this amount of pseuduscalar admixture is 
tolerable at this energy. Whether the same amount of ad- 
mixture is allowable at lower energies and/or is con- 
sistent with the nuclear matter results is still an open 
question at this lime. Up to 500 MeV the effects of N A 
and AA channels in p-nucleus results seem to be minimal 
and the specific choice of relativistic equation used to 
determine the relativistic amplitudes does not appear to 
matter as long as the on -shell NN ampliludes can be 
reproduced with reasonable accuracy. 
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A reblmaic —hiplc w n oit tbcory a t a r — hied in lb c otexi of Tvon gda m c The cfaaic nwhm| — fhmdr far ■ 
famioo proyecutc a infi o a Dine rqnniom Mk m optical poien ual derived boo a fdalivistic odupk toweri n g vertex, k it 
Aow Um the two-body Hrtm an o cmcd vriUi ibe optical pawal i» ike one wid oac pnidc oo in a«« dd io id mer- 
wdnKwcs. 


In almost all the relativistic (Dirac) projectile-nucleus scattering calculations dooe in the past few years, the 
optical potential used has been the relativistic impulse approximation (RlA) ( I ] which is the relativistic ana* 
loguc of the non -relativistic first-order impulse approximation optical potential. The tacit assumption behind 
the use of RlA, together with the Dirac equation, in which the heavy target is taken to be on its mass shell, is 
that there exists a multiple scattering series for the optica! potential (in which the RlA is an approximation to 
the first term). However, neither the existence of the relativistic multiple scattering theory (RMST) nor the 
relation of it to the RlA has been consistently derived. The lack of a RMST not only prevents us from performing 
a systematic study of the higher-order multiple scattering terms, but also from making other corrections, such as 
off-shcll effects and Pauli blocking, in a consistent manner. Therefore it is highly desirable to have a RMST. 

In this work we show that a RMST can be formulated in the context of a relativistic meson exchange model. 
In the following we will consider a scalar "nucleon*' interacting with a spin, isospin zero /I -body target through 
meson exchange. A minimal set of meson exchange diagrams required for any such theory is the set of ladder 
and crossed ladder diagrams. In the limit when the heavy target becomes infinitely massive, this set reduces to 
a one-body equation for the lighter panicle moving in an instantaneous potential produced by the heavier par- 
ticle (the one-body limit [2]). and at high energies gives the eikonal approximation to scattering (3]. We will 
assume that the theory we seek is ere. in which these relativistic ladder and crossed ladders are summed efficiently. 

The contributions to second order in the projectile-meson interaction are box and crossed box diagrams, 
shown in fig. I. Since the target is a complex system with (in general) many closely spaced energy levels and 
continuum states with different combinations of dusters, we assume that alt of these stales can contribute, at 
least in principle, to the intermediate states. These states will be labelled by the index n. with n=0 referring to 
the ground suie. Included in this sum are stales where one, two, and possibly many nucleons are knocked out 
of the target. The elastic scattering matrix for the two diagrams in fig. I is, for spin zero particles, 

1 Preseat address: Department of Physics, Hampton Uaivcniiy. Hampton. VA 23&6I. USA 
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Fig. I. rry— tor— far clastic nutn« of a projcciile 
(mM kw) fioto a ncte (AmUc aolid l«c) is its |ra—d sale. 
Their topas are second order ia Ike woo* (dash hoe) pro- 
jeaie com f Umg g. bur potato to fc ran i ron of the A-bcdj lasgci 
ndeom lo aR order ia g. The todc* ■ fcMd aa arbitrary cx- 
died w of the tar»K. «itb m « 0 toe grotto state. 


Fi*_ 2. Stogalaritirs of ik boa aad crossed boa tor— ia ibc 
cow p fa a 1« plane for a'iJU,+«i. aad f«p"=0. The crtV 
am — 4 ibc meson poles arc to indicate that tbty arc doable polo 
when p-mp'. AR tbc pota a ow aa |t| varirs C as conttog froai 
IhrUTarTarnTflbr nriTlin|i f, tail r , nt am itnr~n 


i/ ? f d 'k^P-p-'P-k^iP-k.P-pWk-p Wp-k) 

,.o * J 

_ ? , f d 'kg*(P-P\ P-q)g^,(P-Q. P-p)J(k -f* )Jjp-k) 

.-o W j (2*) 4 (Pt-*o-i«)lt;(«)-(<V+P, 

where the toul four momentum in the center of mass is P= ( IK 0). q=o’ +p—k and the meson propagators arc 


The meson-projectile vert ides which couple the ground stale lo the nth excited state arc go. and g^o. and 

r 4 w(m I +* 1 )'' 1 , EJk) = (Ml+k } ) t/s . (4) 

where m is the projectile mass and Af„=A/„+4. is the mass or the nth excited state of the A nudeon system with 
excitation energy A., with ^=0. 

The key to our derivation is the analysis or the singularities of the box and the crossed box diagrams in the 
complex kc plane. For this purpose we work at threshold where p=p =0 and H'= \t a + m. The singularities are 
shown in fig. 2. 

(i) If the K, contour for the box diagram is closed in the upper half plane, the pole at A/„ — £„(A ) + m domi- 
nates. This pole corresponds to restricting the target ground slate to its positive energy mass shell. 

(ii) The next largest contributions to the box diagram, and the leading contributions lo the crossed box dia- 
gram. come from the meson poles. Furthermore if the contour for the box diagram and the crossed box 
diagrams are dosed in the upper half plane, these double meson poles at Jt,=m-w nearly cancel. The sum of 
these two terms is proportional to 


I I 2 (Afp — £.) 

Em — Afq+w Ma — E.+w ~~ ut*—{M 0 — E m ) 2 ’ 


(5) 


so that when n = 0 and A/ 0 -*oo. the pole calculated in (i) above gives the exact answer (the negative energy 
poles arc also negligible in this case). If M 0 * oo, an estimate of the ratio of the meson pole contribution to the 
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dominant target pole contribution, (i) above, is fr(m+M 0 )/4inMi. 

(iii) For excited slates of the target nucleus, the pole in the box diagram at M a — E m (k ) + m may overlap the 
double meson pole or even the negative energy nucleon pole at — e*. (This happens when 4=wor 2m. ) Such 
overlaps, which seem to be a manifestation of dissolution singularities 14). introduce spurious singularities into 
the equation. However, all of these spurious singularities are eliminated if the contour is closed in the lower half 
plane. Keeping the pole at k 4t =e L still permits us to separate out the leading term from the box. but in this case 
the cancellations in (ii) will go like (m—e l )/ur. which docs not approach aero as A —oo. However, this loss of 
convergence can be accepted since the contributions from excited states are generally smaller anyway. 

These observations lead us to write the equation for the projectile target r-matrix in the following operator 
form: 

T= l'+ I'GJ T+ VG\, a T , (6) 

where G'J is the propagator for the projectile and the target in its ground state, with the target on its mass shell, 
and (it, a >s the propagator for projectile and excited states of the target, with the pro/eaileoa the mass shell. In 
this equation 

r«jy+r. (7) 

where v ' are the OBE diagrams describing the interaction of the projectile with the ith nucleon in the target and 
I" is the sum of all irreducible terms remaining from the full ladder and crossed ladder sum. Points (i)-(iii) 
imply that V is very small, and if .-f and m both approach infinity the leading OBE terms are exact (2). In 
general cq. (6) sums ladders and crossed ladders exactly if l" is included. 

While cq. (6 ) is an exact formulation of the problem, it is too complicated to be useful. We need a philosophy 
for identifying leading contributions which will be summed exactly and others which will be treated perturba- 
tively. The philosophy wc use is familiar from non-reiativistic multiple scattering theories [5]. The leading 
efiirts are assumed to arise from multiple scattering through the intermediate states in which the projectile 
interacts repeatedly with the same target nucleon. This is important compared to the terms where the projectile 
is interacting with two or more different target panicles since the matrix elements of the latter are proportional 
to (small) correlation functions involving two or more particles. It is not our intention to improve on these 
basic assumptions, but rather to describe how they can be implemented in a relativist ieally covariant manner. 

For clastic scattering, a convenient first step is to introduce an effective potential U (the optical potential). 
In operator form, the r-matrix in terms of the potential U is 

T=U+UGZT, (8) 

where the equation for the optical potential follows from cq. (6): 

1/=I'+KG:, 0 (/=F+1/G;, u K. (9a) 

(/= I + 1 G:, 0 r+ VGZ. n VG mrU V. (9b) 

Eq. (9a) sums all inelastic contributions; (9b) is convenient for projecting the result onto the elastic subspacc 
needed ineq. (8). For larged, eq. (8) is an effective one-particic equation for the projectile moving in a fixed, 
instantaneous field generated by the target. If the projectile has spin J, cq. (8 ) is a Dirac equation. 

To lake into account all the leading effects from rcscaitering from the same nudeon. which controls the strong 
short range NN interaction, wc introduce the multiple scattering scries as discussed above. To this end, separate 
out f from V, introduce a new propagator g, such that 

/'= v> + t'g, r\ (10) 

V’zzV+V’g.V' . (II) 
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ith nucleon in the nuck us, but is otherwise 


(12a) 

where 

(/=!(/'. (12b) 

I 

Eq. ( 1 2 ) is our final result for the optical potential. It gives the exact result for the sum of all ladders and crossed 
ladder diagrams, and in the nonrelalivistic limit (m-»0), 

The propagator g, should be chosen to minimize the contribution from inelastic channels, so that the second 
term on the right hand side of eq. ( 1 2a ) can be treated pcrturbaliv-Iy. If these contributions are dominated by 
one nucleon knockout processes as discussed above and illustrated in fig. 3. choosing g, to be the propagator 
with both the heavy A— I duster and the projectile on mass shell will exactly cancel the dominant inelastic 
contributions from the second term in eq. (12a) and ensure that they are exactly accounted for in the summa- 
tion of eq. ( 10) which produces 1‘. Restricting the A — I duster to its mass shell 11 ensures duster separability of 
the remaining two nudeon system [7], Eq. (10) fori' then reduces, in the NN subspace, to the one particle on- 
shdl (spectator) equation previously introduced by one of us (8], the only change being the shift in the total 
energy of the two-body subspace du** to the motion of the .-I — I duster. 

Finally U can be projected onto the clastic subspacc using eq. (9b). This leads to an effective two-body /- 
matrix for the first term on the right hand side of eq. ( 1 2a ) (denoted by / ) which has both nucleons in the initial 
and final states ofT shdl and is obtained by quadrature from the spectator amplitude t 

i*= p*+ v i g.v'+ v’g.l'g.v* ■ (13) 

Here 0 is the OBE potential with all four legs off shell (unless one of the legs is projected on-shell by g,). In 
applications the first term ineq. (12a) is usually simplified by using the t, approximation and is referred to as 
the Rf A. Our derivation suggests that the full first term of eq. (12a) with g, defined as in fig. 3 is a more precise 
definition. 

This prescription, with some approximations, is precisely the one recently used to calculate p- 4 "Ca scattering 
at 200 MeY (9]. This calculation gives very good results, showing that calculations based on the theory pre- 
sented here can work very well in practice. 

In conclusion, we would like to emphasize that, in the context of the meson exchange model, the projectile- 
nucleus r-matrix docs not readily assume the form convenient for multiple scattering analysis. In order to obtain 
a more manageable kernel and the corresponding /-matrix, we need to consider the explicit cancellations of 


where V'm V /A. Note that g, describes the propagation of the 
unspecified. With these definitions, eq. ( 9a ) becomes 

(/'=/•+ Zr(G% 0 -MU‘+?(i+G:.oZu>y 


" Note Itut i he prescription that the A — I duster be on mass shell has been given in ref. 1 6 J. 
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Fig. 3. Diagrammatic representation of second order mraltcring 
of the projectile with the same target nucleon through interme- 
diate states in which one nudeon is knocked out and the remain- 
ing A - I system is in some cscited slate n,_ Choosing g, as 
shown ensures that these leading terms are exactly included in 
the first term ofcq. (12a), and that the other terms are small. 
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meson poles between the box diagram and the crossed-box diagram. Once the integral equation for the /-matrix 
is obtained, the optical potential can be derived in a straightforward manner. The optical potential can then be 
expressed as a multiple scattering series, eq. (12), and in the impulse approximation the /-matrix associated 
with the optical potential is found to be the one with one particle on its mass shell. 

It is a pleasure to acknowledge helpful conversation with SJ. Wallace, who first alerted us to the problems of 
dissolution singularities. We would also like to thank P.C. Tandy and W. Van Orden for discussions on the 
subject on various occasions. This work was supported in part by the Department of Energy, through CEhAF, 
and by NASA grant NCC1-42 
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Amideuteron-nucleus absorptive cross sections for intermediate to high energies are ca l cu l a t ed us- 
ing an ion-ion optical model. Good agreement with experiment (within IS percent) is obtained in 
this same model for p-nudeus cross sections at laboratory energies up to 13 GeV. We describe « 
technique for estimating antinudecs-nudeus er as sections from NR data and suggest that further 
cosmic ray studies to search for antideuteruns and other antinudei be undertaken. 


The search for antimatter in the form of cosmic ray an- 
tinudei is an intriguing and speculative endeavor. Cer- 
tainly, with present experimental capabilities, producing 
particles heavier than antiprotons in the laboratory is dif- 
ficult at best. Indeed, the question as to whether or not 
one should even search for antinuclei in cosmic rays has 
ton addressed from several different perspectives. One 
point of view, 1 for examp'e, argues (hat theoretical abun- 
dances, estimated from empirical observations, of antinu- 
clei with Z £ 3 are negligible. Although the extreme rari- 
ty of antinudei events may reflect nonconservation of 
haryon number in our universe, the purpose of this paper 
is not to address these issues but, rather, to provide a cal- 
culational procedure for determining whether or not an 
antinucleus has interacted with a nucleus. To do this we 
calculate anrinucleus-nucleus total and absorptive cross 
sections utilizing an optical potential model 2,2 of nucleus- 
nucleus scattering as described below. Numerical results 
for d nucleus are presented to illustrate the predictions. 
Since experimental data for antinuclear nucleus collisions 
are nonexistent, predictions for p-nucleus cross scctior.s 
are made, compared to aval able experimental da* and 
arc found to be in gcod agreement (within 13 percent/. 
Comprehensive tabulations of the predicted antinucleu.- 
nudeus cross sections are published elsewhere. 4 

For the scattering of composite nuclei, a general 
multiple-scattering theory (neglecting three-body interac- 
tions) has been developed by Wilson.* The series reduces 
to the usual Watson form when the projectile is elementa- 
ry. Through the use of the impulse and closure approxi- 
mations, a simple, folded, optical model potential wjs de- 
rived 2 as 

Wz)=A,A r f dhpriz) 

X f d'yppix +v+z)7le,y) , (1) 

where j is the NN kinetic energy in the c.m. frame, y is 
the NN relative separation, p T and p f are the target and 
projectile number density distributions normalized to uni- 


ty, Tfcjd is the energy-dependent constituent-averaged 
two-nucleon transition amplitude obtained from scattering 
experiments, and Ay and A T are the projectile and target 
atomic numbers, respectivdy. With no renormalization jr 
parameter adjustments, this optical potential has been pre- 
viously used in a Wentzd-Kramers-Brillouin (WKB) for- 
mulation to obtain excdlent agreement with experimental 
elastic scattering differential, reaction, and total cress sec- 
tion data at energies lower than 23 MeV/nuckon.* More 
often, however, this optical potential approximation is 
used within the context of an cikonal formalism to predict 
nucleon -nucleus, deuteron-nudeus, and nucleus absorb rive 
(inelastic) cross sections to within 3% for energies higher 
than 80 MeV/partide and to within 10% for lower ener- 
gies. 2 * 2,7 From eikoital scattering theory, the absorption 
(reaction) cross section is 

2>r J o " |l-expl-2IinYU»)l| , 'db . (2) 

where the complex phase function is (with Ii= 1 ) 

AT(*)=— (2 kr'f m V(b,z)dz , (3) 

* -• 

with k the projectile momentum wave number and b 
denoting the impact parameter. The reduced potential is 
then obtained from the optical potential as 

lHx)=2mAyA T iAy+A T )~ l W(x) , (4) 

where m is the nucleon mass. With : n the eikona! context, 
this model is similar to the comparable, but alternative, 
Glauber iheory formalism which has beer, extensively 
developed by Franco and collaborators.* Aside from the 
improved convergence to the exact multiple-scattering 
series by the Wilson approximation* (due to difTer<*nces in 
higher order terms), the Wilson propagator 2 also includes 
target recoil and terms to ordet kL 
In order to apply Eqs. (1)— {4) to antinuclcus-nudcus 
collisions, several assumptions, other than the applicabili- 
ty of the underlying composite-particle multiple-scattering 
formalism, are necessary. First, we assume that the num- 
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her density distribution for an antinuclens is i d entical to 
that of its "normaF nudess co unt erpart. Hence, aside 
from the overall sign of the charge distribution, the an- 
tinucki charge d i stribution parameters and functional 
forms are assumed to be identical to those of the corre- 
sponding nudesr species. We then extract the numb e r 
densities p A from the corresponding charge densities pc 
obtained from dectran scattering experiments, by assum- 
ing that 

pt<r)= /, *<'■>„ rWV . (5) 

where the mtdeon charge distributkm is t ak en to be the 


p^r)=(?/2»rii , ' 7 expC-3r 1 /2ri l |. (6) 


with a nackor rw mean-sqnare charge radius set equal 
to the protou value* or 0 l* 7 fnt. 

For elementary proj ecti les . subsutoitou of ( 6 ) into (5) 
yields a deha function far p A . For nu lei or ant i n udci 
with A <20 we use aa harmonic osc ill a to r form fb^ts 


pc<r)^pt 




expi-rVa 1 ) . 


(7) 


where the charge distribution pa ramet ers j and a are bat- 
ed in Table I. Inserting ( 6 ) and (7) into (5) yields 


^{r)=(^ , /«r , ICl+C3r/2)-C3r« 2 /*s I l 

+(r«V/16s < )lexp(-r I /4» 1 ). (« 


where 

s , =(« 2 /4)-(rJ ( /6) . C9» 

For aedri or antinudri with A £ 20 we choos e a WootH- 
Saxoa form 


TABLE L Nudesr charge dnt ri butror p a ra m eters from dec- 


traa scancriag data. 


Kmfcm 

Dotribtioa* 

r ihw) 
or t (fmXWS) 

a (HW) 
ar it (fmXWS) 

*H 

HW 

0 

1.71 

*He 

HW 

0 

1.33 

’Li 

HW 

0l327 

1.77 

*Be 

HW 

0.611 

1.791 


HW 

axil 

1.69 

,J e 

HW 

1-247 

1.649 

M N 

HW 

1.291 

1.729 

“O . 

HW 

1.544 

1.133 

"Ne 

WS 

2.SI7 

2.74 

"A! 

WS 

2.504 

3j05 

"Ar 

WS 

2.693 

3.47 

**Fe 

WS 

2.611 

3.971 

**Cu 

WS 

2.504 

4.20 

"Be 

WS 

2.306 

4.a04 

•"A* 

WS 

2.354 

5.139 

,M 3a 

WS 

2.621 

5.61 S 

*"Pb 

WS 

2.416 

6.624 


'The lumccic wdl (h .Lstributioa is used for A <20 rod 


the Wood»-Sa on (WS) dnlnbutioo for A ? 20. 


Pc=Pol 1 +exp((r-/t )/c]| 1 , (1« 

where Jt is the half-density radius, and the surface dif- 
fuseness c is related to the charge skin th i ckness 4 
through 

t r = 4 . 4 ' . ( 11 ) 

Values for it and t r are also listed in Tabic I. Inse r ting 
( 6 ) and <10* into (5) yields, after some simplification . 1 a 
number density p A that is of the Woods-Saaon form with 
the same it, but different overall normalization p% and 
surface thickness. The latter is given (in fm) by 

^J=3.OCr N iln((30-l)X[3-/»]^ , , (12) 


0=tM1.5Ar H /t c } . (13) 

In all cases the densities are normalized to unity. 

For the aat in n d eo n au de ou (NN) transition amplitude. 
7. ire assume a farm 


JlejF)= 


1/1 

~ rJrl(alf)+iK2rWr)l" M 


XexpC-/*/2B(e)l. (14) 

which is obtained by taking the Fourier transfonn_of the 
usual unclear amplitude. 1 * In Eq. (14) <r M is the NN total 
cross section, crie) is the renl-to-imaginary ratio of the for- 
ward amplitude, and Ble) is the NN slope parameter. 
With this choice for the *wobody amplitude we effective- 
ly treat antiun.~leas nucleus scattering as a purely comple- 
mentary system to nucleus nucleus scattering. Values for 
r M ate param etrized in terms at the incident pa r t icle 
■ om entu m Pby 

o — =(61.2 mb) +(53.4 mbGeV/e)/P 

— P(mbGeV/c) . (15) 

which is c modification at the expressior given in Ash- 
ford rr of . 11 to extend the cross sections to IS GeV. 
Values for the slope parameter, displayed in Fig. 1. were 
taken from the Paris NN potential 12 and from Block and 
Cahn . 11 Values of the parameter a are nor presented for 
these results since only the imaginary pan of the transi- 
tion amplitude is used to calculate absorption cross sec- 
’iora. Details of the values used for the total cross section 
estimates are found in Ref. 4. 

Predictions for p-nuckus and d-nuclcus absorption 
cross sections are displayed in Figs. 2—4. Also plotted in 
Figs. 2 and 3, for comparison, are data from various ex- 
perimental measurements . 14-1 4 Clearly the agreement be- 
tween theory and experiment is good since the maximum 
cross section difference for any of these results is less than 
15 percent. Typical differences between theory and exper- 
iment are 5 percent. Figure 4 displays predicted d-nuclcus 
absorptive cross sections for the range from 50 
MrV/nudcon to 15 GeV/nudeon. These are provided in 
the event that techniques for producing antideutcrons in 
the laboratory may become available in the future. Simi- 
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may also contribute to tbe usefulness of this collision 
mood. Indeed a primary cosmic antitriton event was re- 
ported a*, the 18th International Cosmic Ray Confer- 


up to Fe-Pb are tabulated in Ref. 4. For all collision 
pairs, tbe cro ss se c ti on s , as in Fig- 4, are smooth c ur v es 
displaying no complicated structure. For optical roodd 
calculations, however, this is not entirely unexpected since 
it merely reflects the averaging nature of the calculation 
and tbe smooth fractional d epen d en cies of tLe input NN 


FIG. 2. Theoretical p- nucleus absorption cross sections as a 
function of incident kinetic energy. Experimental values were 
obtained from Refs. 14—16. 


FIG. 3. Theoretical p-Pb absorption crass sections as a func- 
tion of incident kinetic energy. Experimental values arc taken 
from Refs. IS and 16. 
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FIG. A Theoretical predictions for d-Mdm abaorpervr crass 
sections as • function of mridett kinetic onp. 


We have abo per formed an analysts to investigate the 
iwiiiiiUji of oar predictions to the magnitude of the dope 
parameter. This was doae by recalculating aatiproton- 
aad Ft arnicas crass scctioas for the cases where the NN 
slope parameters are tritea ia a rather md hoc fashion as 
2JK*) awl *(e)/I These slope paraaieter ranges are 
displ a yed ia Fig. 3 along with the Brookhaecn measure- 
ments." As anticipated, are found that cotlisiaas involv- 
ing light Badri me far more sensitive to chang e s m the 
slope parameter than collision s involving heavy nuclei, 
especially at lower energies. For instance, with an as- 
sumed slope parameter of 25, we Had that at 100 
MeV/andeon. increases by 21% for p-C (629 mb vs $20 
mb},_bat oa!y by 10% for p-Pb (2528 mb vs 2303 mb). 
For Fe-Pb at 100 McV/nuckoa. the increase is 6% (5886 
mb vs 5542 raid. At IS Gcv/aadeoa, the i ncreas es for 
these same collision pain are 13%. 7%, and 4%, respec- 
tively. If the slope parameter is halved (to 0.5 Bi. we Had 
that a* for the isme three collision pain decreases by 
13%. 6%, Md 4% at 100 MeV/audcon. and by 9%. 4%. 
and 3% at IS GeV/nndcon. Clearly these absorption 
crass sections are not very sensitive to large changes in the 
slope parameter. 

la summary, we have employed a simple optical model 1 
from nodens-auckus scattering theory to intermediate en- 
ergy antinaricus-nadeas collisions. The only new inputs 
required to complete the calculations are the experimental 
NN elastic scattering parameters (total crass section, elas- 
tic slope parameters, and real -to- imaginary ratio of the 
forward scattering amplitudes). For the energy range con- 
sidered herein tiOQ McV/nodeon to 15 GcV/nudeon), the 
akonal formalism is certainly adequate. Although these 
methods could be extended to even lower energies. 7 other 
more suitable methods'* air currently being implemented 



t aw awn 


FIG. 5. Range of NN dope parameters need ia the crass sec- 
tion mainty studies. The experimental vaiaes are taken from 
Bnukbra acagwaaut SRef. 1 1). 

elsewhere. We have included charge exchaagejmly as it b 
included in the determination of the Paris NN slope pa- 
rameters. We have not considered pioa propagation ef- 
fects in the target nucleus, two-body correlation functions, 
and have ignored possibly important ^pin and isosptn ex- 
citations of the target. 1 * Nevertheless, we find that the 
model gives rather good agreement with available 
p-eudem dataju intermediate energies and expect that 
the predicted d-nedcus absorption cross sect ions are 
reasonably accurate (certainly within IS percent). While 
are realise that the product ion of enough antideutcrons in 
the laboratory (LEAR for example) to produce a beam is 
in the distant future, the production of such a beam could 
open up new areas of research since, for example, high 
temperatures in nuclear matter, may be achievable 20 with 
antideutcrons of momentum greater than 2 GeV/c. For 
now. cosmic ray studies appear to offer the best 
for detecting and studying d and antinuclei interactions 
such as the reported xntitriton event- 17 

The authors wish to thank Carl Dover, Dan Strotunan, 
Robert Welsh, and Rolf Winter for useful discussions. 
Support from NASA under Grant Nos. NAG 1-477 
(W.W.B.) and NCC 1-42 U.W.NJ is gratefully acknowl- 
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Abstract 


Resu t are presented of an analytic relativistic calculation of a OBE nucleon-nucleon (NN) 
interaction employing the Gross equation. The calculation consists of a non-relatmstk reduction 
that keeps the negative energy states. The result is compared to purely non-relathrbtic OBEP 
results and the relativistic effects are separated out. One finds that the resulting relativistic 
effects are expressable as a power series in r t • rj that agrees, qualitatively, with NN scattering. 
Upon G- parity transforming this NN potential, one obtains, qualitatively, a short range FfN 
spectroscopy in which the S-states are the lowest states. 



In the quest to understand short range nuclear forces, many nuclear theorists have embraced 
QCD and have made progress toward extracting the roles that quarks play. This approach 
is further fed by the observation that traditional non-relativistic meson exchange models have 
difficulty providing any new information on these short range forces. In this paper, new results for 
the short range interactions are presented that arise from a relativistic view of the N-N interaction 
that contains negative energy states. Preliminary and partial results were presented earlier^ 1 ) 
This relativistic approach is not new; it was introduced by Gross in 1974.W In this present work, 
that pioneering work is revisited and expanded upon in such a way as to provide further insight 
into the nature of the short range forces harboured in the relativistic wave equation. It will be 
shown that, for the One Boson Exchange Potential considered, the short range contribution can 
be expressed as a power series in - rj, the Nucleon-Nucleon isospin operator. This short range 
contribution is interpreted as a relativistic effect and is a direct result of coupling to negative 
energy states. It will further be shown that, as a consequence of this relativistic effect, the G~ 
parity transformation of the elastic NN potential^) gives rise to a new level ordering prediction, at 
short ranges, for the Antinucleon-Nucleon (N N) interaction. The final result of the analytic work 
presented below is a qualitative description of the (NN and NN) interactions, as no numerical 
values for the exchanged mesons’ masses or coupling constants are employed. A word of caution 
is needed, however; the NN interaction as presented in this paper is not antisymmetrized. The 
impact that this may have on the short range NN contributions presented here is uncertain; 
however, it is assumed that antisymmetrization does not apply to NN. 

The starting po : nt is the Gross equation^ 2 ) for the NN system written as: 

(few?) — f (i) 


where 


total 4 — momentum = P = (IV, 6); W = total rest mass 

(M + ^ + 


G^fy^k.W) = 


2E k W{2E k - W) 


k = {k 0 ,k)\p = (po.p) 

w w 

ko — E k — ; po = Ep y 

E\ = Af 2 + k 2 ,Af — nucleon mass 


( 2 ) 


f is the covariant two body vertex function that is state dependent, and C is the charge 
conjugation matrix. To facilitate making a non-relativistic reduction in order to expose the 
analytic structure of the interaction, it is useful to write 


1 



( 3 ) 



By introducing 


[2W (2jt) s J j 
-A/ 

Vr,W (2W(2*)»|§ 


4 r) (p)«l ~ <} (-p)rc)^(p) 
^(2£, - »T) 


4 r W~ <} (p)(rcW(p) 

E f W 


(4a) 

(46) 


where V> + and V ore the positive and negative energy momentum space wave functions 
respectively, the following set of coupled equations can be extracted: 


~ W *rM = -/ (0 )} <«> 

The V' ++ ,V +_ , V~ + t and V are related to the one particle on the man shell interaction 
kernels V^ vt p by: 


v++ = )4‘ , (-f)v w .^.M(t )«!;''(-* ) (7> 

v + - = )«£■*(-? )V.^“ ,r) '< k ' >•!?(* ) w 

v ' + “ )«£’(-*) w 

v ) V rS,>-s u s (i )»*7‘ ’(* ) (10) 


Of course, the non-antiaymmetrized V^i represent meson exchanges and, as is customary, these 

interactions will be approximated by single boson exchanges; namely, x, a, p, and u. One notes 
that there is no concern at present for the numerical values of the masses and c upling constants 
of these bosons. Thus there is no concern that the interaction not reproduce the NN phases, 
effective ranges, etc. One quite simply wants to compare qualitative features of the relativistic 
interaction to that of the non-relativistic interaction. This is performed by, essentially, subtracting 
the non-rel&tivistic interaction from the relativistic interaction presented here. That is, the limit 
as r — ♦ 0 is taken. What remains from this procedure is what one considers the relativistic effect 
or simply, the interaction difference. To arrive at results that can be treated analytically, a non- 
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relativistic reduction is performed that keeps the negative energy states. Having stated this, one 
continues with the calculation. 

Employing expansion approximations such as *** * 1 « 4 

with = 4^l(l~ 4^71 — 4^1 + )* After quite a bit of algebra, equation 

5 is reduced to 


x { 1 - [*•»*• (2ft* ) + 2* ’] *£, (* } + jjy • (rf ■ - f )+,/ (* ) } (11) 

for the sigma exchange only. Equation 11 is the result of keeping the lowest order °f * or -j^- 
compared to the leading term. Equation 11 as well as the other boson exchange contributions 
can now be Fourier Transformed to configuration space. (Similarly treated is equation 6.) The 
motivation for going to configuration space is the ease in which non-reiativistic and relativistic 
contributions can be compared. Traditionally, non-relativistic potentials are always presented in 
position space. Keeping in mind that we seek only qualitative comparisons, we then transform 
our momentum space reduction into position space. 

The resulting configuration space coupled equations are written: 


-(^ + ^+e)* + (?)=-»' ++ (?)« + (r-)-K + -(f)*-(r) (12) 

-2 = -V-+(f )#+(r)-r--(r )*-(r ) (13) 

The results of Equation 12 and 13 are not new. Gross presented these equations^) without the 
quartic derivative operator. 

The potentials V ++ , F +_ , V and V are: 


V++ = «e + o\ • <rju„ + S 12 ut + L • S*LS 

= !fLli + vr + llM£l Vf+iei.Vv' 


v-+ = (v+-) f 

, L • 5 


= V. 


* Mr 


V* + V a w - 


L»SV" 


+ 


Ol • <72 V/ + SilVf + V> - 


(i + JA >) s 


Mr 


L*S V[ 


+ higher order terms 


(i<) 

(15) 

(16) 


(») 
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where 


and 


u e = Vf + V" + f r 2 V' 

■ •» + M J v/] + 1 + W’C 

«r = n • i [*7 - (i + «,) ! ] - (i + M**? 


«£S = ~ 


m, 

Afx 


P7 + (1.5 + 2 AO Vf + t! • rj((l.5 + 2*>) V/] 


(18) 


* = 

*7 = 

*r = 
*r = 


KTM 

*7M 

V7(r) 

V/(r) 

V/(r) 

V/- W (r) 

vr W 

V/' W (r) 


rm r 

~ r i • ^n r 

-(i + *.)iT-n*»»(i + *p)v/ 
^[vZ + C + n.^v/j 

—a* 


m * 

17 m ' “7 


where a = 

m r 


& 

ml 

^ I 

Ax 

2 M 

X ' 

9± 

m\ 

e~ x 

Ax 

12M 2 

X 

Xgl 

Ax 

m l I 
2 M ' 

K) 

9± 

ml 

( I+ r 

Ax 

12M 7 

9p,<*> 

m r 

e -/>*,W* 

Ax 

z 

9p,u 

ml 

(i 

Ax 

M 

* V 

9 1 

_3 

e-A i 

Ax 

12M 7 

z \ 


R) 


trim mu , 

p = — ; w = — 

m r m x 


(10) 


These relations are well known^ 2 ) and one can verify them. The next step is to uncouple Equations 
12 and 13 to obtain a single Schroedinger like equation. One finds that 


" (if + *) * +w = “ ( v++ - v *) * + w 


( 20 ) 
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where 



y - 1 y + ~ * y-+ 

4M? 

(21) 

This is the non-relativiatic reduction to be examined. In all of the work that follows, V has 

been neglected. Thus 

V,- v< + l v+ 'l’ 
R 4Af s 2M 

. (22) 

One can show, through the employment of spin and angular momentum “aerobics” 

that equation 

20 is equivalent to 


(23) 

where 

/ 2\ 



-('*?) 

(24) 

and 

= K + <?i • o 2 V M + 5 i 2 V$ 12 + L • SV^5 + £• DV^q 


V T 

(25) 


The V’s are the same as those described in reference 2 and are a convoluted arrangement of a, p, 
a, and w potentials. 

For the next phase of the calculation, one can proceed either from Equation 25 or from the 
potential found in Equation 21. Proceeding from the former choice, one finds after perf ming 
some algebra and keeping only the largest contributions as r — ♦ 0: 

V c = a c + b e r t • r 2 + c c (ri • r 2 ) 2 - a' - 6'rj • r 2 - c' fa • r 2 ) 2 
V LS ~ ~ a LS ~ t>LS T l * r 2 ~ c Ls[ T l • T 2 ] 2 
Vsi2 = ~ a 512 ~ &S12 T 1 • t 2 + c 512( t 1 • T 2)* 

= a„ + b mTj • r 2 + c„Ti • r 2 ) 2 

V LD = a LD + biD T 1 • r 2 + c LD{ T i • r 2)* ( 26 ) 


where the a’s, b’s and c’s are positive definite and, to the leading term, are- given by 


«« = lFV.“' 

ot,=GV? 

= iK,”’ 

“LD = ifv;v? 
*S1J = f V“ J 

= W? 


l (l2+31 > ) v p vw 

0e — jJ/r “o *o 

= 3GV'V?' 
b„ = £{1 + 3 k,)V'V? 
h D = FV“ V ! 

*su « t(1 + 3k,)V'V“ 

ii5 = 3FV“V/ 


Ce — (4 + 3fc, + 
4 = 3CV' , V r " 

e.. = f V/* 

= j''/ 

c s n = n( 3 + 
C I5 = f y/ 


and 


(27) 
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P — 


— * Cl 


A/ 3 x* 

H 

V 

= 4 

4 w 

e~ Bx 
m r 

X 

; B- 

= a,p,u 

m a 

m, 

; p = — 

; w 

mu 

m r 

m r 


m r 


= m»«y 

4 mj 


£>r = i + 


( 28 ) 


One notes thatjin this limit the pion contributions can be neglected compared to the other terms. 
To obtain the NN potentials for small distances, one G -parity transforms the NN potentials of 
Equation 26. This effectively changes the sign of the omega coupling constan' and, thus, changes 
the sign of the corresponding coefficients. 

The final ingredients that we need before making concluding remarks are the spin matrix 
elements; all but L • D can be found elsewhere^*) and the L • D matrix elements are found in 
reference 2. For 33 P 0 NN, one finds L — —2, Sjj = —4, a\ • <tj = 1, and q • rj = 1. For 
13 fo N N, one finds that only rj • rj changes; rj • ry= -3. Both of these NN*P 0 states have 
the same Z* I). Making the substitutions into Equation 25 gives the qualitative resultr ' he 
I3 Po potential lies higher than the 33 P a ; a result in agreement with the numerical work ■ ■ •' 

researchers.^ 4 ). Furthermore, through similar arguments, one finds that the 11 S’« NN ies ?- 
than the 13 P 0 NN; an tmexpected result. Finally, it is clear that all isoscalar NN pot e 

more repulsive than their isovector counterparts. Hence, 11 5 0 > J1 S„ P, >» P„“ D, , 1 D,, 
etc. 

These qualitative results should be unaffected by a more complete interaction since it is well 
known that the omega meson exchange dominates the short range interaction. The omega meson 
exchange is included explicitly here. More complete interaction models should vary only in their 
quantitative results such as the amount of energy level shift. It is not clear if the results presented 
here will affect the pp Coulombic states widths. Although, theoretical aDproaches generally ‘cut- 
off* the pp interaction inside 1 fm, investigating how il ® relativistic effects affec Coulombic state 
widths is worth pursuing. 

In conclusion, to obtain analytic results, the Gross equation was examined in a non-relativistic 
reduction of the NN interaction that keeps the negative energy states. The NN interaction was 
chosen to be a one boson exchange consisting o r x, a, p, and u>. The reduction was then applied to 
the real part of the NN interaction via G-parity. To get a qualitative xeeling for what coupling to 
the negative energy slates provides, a short distance limit was taken. One might expect that the 
difference between relativistic and non-relativistic theoretical descriptions would_ show up at short 
distance. This work finds that indeed that is the case; for the 11 S’© has a real NN potential that 
is more attractive than that of the 13 P 0 ; a result rather different from reference 3. The fact that 
this is the case at very short distance for this work or any other work may be worrysome since 
annihilation was not taken into account. On the other hnnd, there is no conclusive evidence that 
annihilation contributes any more than giving the states widths. Furthermore, this “relativistic 
effect” may start to be evident at ranges as long as 0.4 fm in some channels. An effort is already 
underway to include annihilation in order to calculate cross sections and other effects. One final 
note is that the level orderings are directly related to the isospin coherences of Equations 26 and 
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from a purely non-relatiristic viewpoint this can be thought of as a result, in part, of adding the 
contribution of a Z graph. 5 

I wish to thank Franz Gross, Carl Dover, J.-M. Richard, and Rolf Winter for valuable 
conversations and CEBAF for its continuing hospitality. This work is supported in part by 
.ional Science Foundation grant RQ-8704038 and NASA grant NAG- 1-447. 
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Lattice piigc calculations ( I) for sialic (heavy) quarts sup- 
port the notion that the urterquark potential in quantum chrtv 
modynaniics (QCD) behaves as V (r) — Ar for lap: r. Indeed, 
the linear potential has loop been used m phenomenological 
.mnrdativislic quark models of baryons and mesons (2. )|. 
Meson spectroscopy in particular is successfully described by 
a linear potential at large r. modified by spin- and ctdour- 
depende— Coulomb forces — small r. Mast calculations >-rth 
the linear potential ate carried out in coordinate space. This is 
the simplest procedure for heavy -quark systems, which can per- 
haps be considered as nomelativisiic. however far light-quark 
systems it would he desirable to have a relativistic treatment 
Bound-stale equations in relativistic systems (4) arc penerJIy 
much easier to solve in momentum space, and thus we arc led 
to consider, as a starling point for the relativism ease, the 
Schrodinger equation fur two scalar particles interacting hy a 
linear potential. The methods developed will generalize rela- 
tively straightforwardly to irlalivistic treatments. 

To summarize: here, we treat the Sdwfidingcr equation for 
a linear r-space potential. The method is for die most part 
straightforward, the only difficulty arising from the singularity 
of the kernel — the origin of momc—um space. Previous treat- 
ments t5) have UMcniy been appro. --male in the ser.vc that the 
singularity was handled hy screening ;ne r-space ponlcntial: 

|l| V(r) -kre* 

What has perhaps not been generally appreciated is that the 
limit t) — ► 0 can he taken analytically. Previous treatments keep 
the parameter n finite, leading to some uncertainly as to the 
nature of the calculated eigenvalues and wave functions. In this 
connection, recall that the screened linear potential docs not. 
strictly speaking, possess true bound states, instead it has scat- 
tering resonances, which for low energy approximate the bound 


stales o< the unscreened potential We vrill extract the limit of 
zero s c ree ning analytically, using a subtraction technique. The 
resulting subtracted integral equation is relatively easy 10 han- 
dle numerically. An alternative procedure, nrt c.np fo y w g any 
subtraction, and leading to a different integrodi ITctent lal equa- 
tion is presented in ref. 6. Our appro a ch is easy to i mple m ent 
and generalizes without difficulty to higher partial waves. The 
Schrodinger equation far the *lh partial wave is (with the 
inhomogeneous term aiteady omitted, as it will not contribute 
to the bound stales in the limit of zero screening) 

121 ^4*/p) + J V/pj>’)+fp’)p' 3 dp' = £+//>) 

Here |i - m t mj(m, + or.) is the reduced mass and V,. given 
by 
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is the Ith partial-wave component of the Fourier transform of 
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The variable v is given by: 
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Q'jL?) and (T,l v) arc (he first and second derivatives (with 
respect to v) of the Legendre function of the second kind. To 
illustrate the method wc specialize to j- waves. where we find 
hy contour integration 

M J vjf/vV’V = £ f «v»' [(?^r> * 

Note that when q = 0. @jfr) and ^Jr) have dnuhtc and 
rpadruptc poles, respectively, at p' — p. so that their integrals 
do not csrst separately. Nevertheless, the two terms added 
together produce a function with an integrate singularity. This 
is illustrated in fig. I . which shows the kernel as a function of 
p' lor fixed p. One observes that there is a central maximum 
at j. - p with height scaling as I/q 1 . flanked bjr two minima 
at p" ~p i 2q whrac heights also scale with 1/q*. The integral 
vanishes |6| and dus allows us to rewrite the Schtudingcr equa- 
tion in subtracted form 


m !f-+jr) + 

'■a 


; / f 


Q'Jy) (S ri 

ipp'f ^ </*»'>' 


x I +jp') - dp’ = 


The limit q — • 0 now exists, and may he extracted hy splitting 
the region of integration to isolate the singularity. We write 


( 8 ! 


j dp' i+j/’ t - +Jj»l 

. r * p * r 

Ja J r 4^ 4. i 4^ 


= a + « + c 

The limits /» * 4q arc chosen so that all three extrema of the 
kernel lie in the middle region B. The explicit forms of the 
Legendre functions arc 


COn = 


I 



Hr. I . The s ingul arity str u ct ure of the kernel is 
i) = 0.075 with fixed p = 2. 


for fi 


** [</»’- /») ; + v + </»' + py +t] 


and 


(?3 .n = + p’ : + q : > 

PP 


[</»' - py * t + f p * + d : + v ] 


It is clear that for p' * p. as is the case in the integrals A and 
C . the limit q — » 0 is innocuous, and may be taken immediately, 
indeed one has 




lim |it + C| 
u -*> 


where P denotes as usual, the Cauchy principal value of the 
integral . which has been made well defined by the subtraction. 
The term B must he handled with care, however, since p' = p 
inside the region of integration. Assuming 4*n') is analytic in 
ilk- RcighhortkKk! of p. and nuking an obvious change of 
variable we find 


1 101 


8 • 'z r , 4 '*[? rv 4 ♦%=][*' 4 S 4 ‘ 4 ]] 

+ I q‘|(r + pV + p 1 + q ? J| — , + ; — ~r ; | fc«|»' + ~ +’+... j 1 1 = Inn /II + Inn B2 

L L» t q (» t -/»» 1 tv J l 2 J J J 


Scaling out 4q then results in 
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N 



iiu.1 

8 

10 

12 

14 

16 

18 

F. 

5.17.1 

5.172 

5172 

5.172 

5.172 

5.972 

5.971 

E, 

10 468 

10.444 

10.443 

10.443 

10.44.3 

10.443 

10.441 

E. 

14.381 

14 114 

14.111 

14.104 

14.104 

14.104 

14.101 

F. 

IX. Mb 

17.452 

17 378 

17.341 

17.335 

17.335 

17.335 

E % 

11.402 

21 125 

20 *47 

20.351 

31.214 

20.293 

20.291 

F. 

27.20ii 

25 6* » 

23.440 

2».2*t 

234172 

2 3 U5t 

21 046 

F, 

33.032 

3I.2WI 

27 274 

26.051 

25.842 

25.648 

25.646 

F . 

44.374 

3*1.224 

32.113 

21.032 

21.789 

27147 

28 111 

E, 


40.311 

38.146 

33.051 

31.561 

30.194 

30.488 

F m 


51.774 

45.304 

38.0b7 

34.428 

33.340 

32.761 

F„ 



44 440 

44.286 

38.517 

36.489 

34.172 

F, ; 



58.588 

51.813 

43.615 

37.301 

37.101 


The contribution of the second term ui HI ctcariy vanishes ante itisM singular ai p = /». the analysis of HI is similar, and we 
conclude dal H lends hi zero. There I ore the limiting firm of die conation is 

I * —I f +Ap) *>] r *r P v IW> - WpH = 

ip. r »/»- J- \ip--pr J 

We now discuss die numerical volution of |I2J. winch is nol yet a completely trivial maner. smee carr must he taken 10 obtain 
die Cauchy principal value. In this respect there is a difference between the linear potential and the Coulomb potential, the latter 
giv.ng rise to a logarithmic singularity. Is* the Coulomb potential, the method used in the lucrative (7) is to write the Coulomb 
analog of 1 12] directly, fir c cample, using Gaussian quadrature. as a matrix equation. Since the singularity is only k^garithmic 
diis method is successful fir die Coulomb potertial Ik-re. such an approach is mu feasible. Instead, we expand in a suitable 
set of basis functions 

X 

1131 Wp) = X 0.<P> 


Inserting this expan Jon in 1 121. muhiplying hy p : ajp) and integrating over p. we obtain 

ll-M £ C.jJ ~xjp)gjp)fip + ~ j [ ^ ] CjpH.C.fp'l “ JT.f P>l dp’ «V» = * £ C. J p' x j p)sj p) dp 


which is just the matrix equation 

1151 £ = E V o__C. 


The double integral over p and p* is performed by changing to 
variables ( p' + p) and ( p' — p). The singularity is in the 
integral over ( p' - p). so this is carried out fust using Gaussian 
quadrature with an even number of points. This type of inte- 
gration yields the Cauchy principal value automatically (X). A 
convenient set of luncliuns j( /») is 


|I6| K.(pl = :- 


(i»Vm* + p 4 


where N is the maximum number of functions used in expansion 
| I V. Figure 2 is a 3/3 plot of the kernel oi l 14|. showing clearly 
the cancellation that leads to the principal value. Using the 
above method, we have calculated both eigenvalues and eigen- 
vectors. In Table I the first 12 eigenvalues arc listed. We used 
m, = m, = 1.5 GcV and the string tension X = 5 GeV’. One 
can see that the lower eigenvalues converge nicely as the num- 
ber of functions is fie teased. We compare these with the eigen- 
values obtained from a coordinate space calculation (integrating 
the equation out from r - Cl and in from large r, and matching 



Fa;. 2. A ihrre -dimensional figure if the subtracted, regulated inte- 
grand: q = 0.075. The cancellation that produces ihe Cauchy prin- 
ciple value is evident. 


the logarithmic derivatives at the . issical turning point), in 
Table I. The calculated eigenfunctions also agree with thc 
coordituie space calculation. 
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In cuocluMua. we have Imlcil the problem of two nonrela- 
livink. scalar panicles interacting via a linear potential in 
momentum space. The relevant Schnidingcr equation has a sin- 
gular kennel. We have shown how after regulating the singu- 
larity by exponentially screening the r-spocc potential, the 
seventy of the singularity can he induced by a suitable sub- 
traction. and the limit of «m screening extracted analytically. 
To the best of out knowledge, this point has not been generally 
understood in the literature The limiting form of the equation 
has been treated numerically, and the results are in good agree- 
ment with more straightforward coordinate space calculations. 
Relativistic equations involving linear potentials involve sim- 
ilar singularities, so that the methods de ve l o pe d here will he 
a pplicable. We intend to study the relativistic quark-anliquari 
problem in the future. The method presented here can be 
generalized to higher partial waves without undue difTicuhy. 

Acfc nwwfadgcmcnts 

We are extremely gratctul to Ran/ Gross for his generous 
contribution of essential ideas during all stages of this work. 
We would also like hi thank Warren Buck and J.W.N. would 
like to thank Rank Cucmutla and Barry Ganapoi for useful 


conversations. K.M.M. and J.W.N. would like to thank the 
Continuous electron Beam Aceckrahir Facility for its hospi- 
tality. This work was supported m part by the Depart mem of 
Energy ihrough Continuous Electron Beam Acce l er a tor Facility 
(D.E.K.). and hy NASA giants NAG I 11.14 (J.W.N.) and 
NAG- 1 -477 (K.M.M. ). 

1. I IV Stack. Phys Rev IV. Part Refch. 27. 412 IPW3V. A 
liascufnu* aud P. IlncUiar. Am. Rev. Nucl. Pul. Set. JJ. 5.W 
I PRO). 

2. li. hktuen. K. &«(iial. T. Kinashiu. K. D. Lane, and T. M. 
Yan. Phyv Rev. IV Put. KcUs. 17. WHO | IS7K); 21 . 203 1 ISHD) 

3. S. Godfrey and N. hpr Pbys. Rev IV Pan. Helds. 32. IIW 
UWl 

4 W. W. Back and F. Gross. phys. Rev. C NuH Phys. It. 23ftl 
tlHT9);F. Gross. Phys. P.cv. ISS. 144* tlUWI; II. Crater and P. 
Van Ahlinc. Pbys. Rev. IV Pul Hchh 37. PK2 tFMni. 

S. A. B. Hcn ri q u cs . B. H. Kctlcfl. and R. G. Mnatuec. Phys. 
Ixtl 44R. K5(IH76) 

6 J. K. Spence and J. P. Vary. Phys. Rev. O. Pan. Fields 35. 2191 
(PRO) 

7. Y. R. Kwm and F. Tahakia Pbys Rev. C: Nad Pbvs. II. U32 
|P»7K); R. II. landau 27. 2IHI |P*3». 

S. I. It. Sloan. J. Cumput. Pbys. 3. 332 1 IUM). 



